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Introduction 


Pipe, What is It? 


Without a doubt, one of the most efficient and natural simple machines has to be the pipe. By 
definition it is a hollow cylinder of metal, wood, or other material, used for the conveyance of 
water, gas, steam, petroleum, and so forth. The pipe, as a conduit and means to transfer mass from 
point to point, was not invented, it evolved; the standard circular cross sectional geometry is 
exhibited even in blood vessels. 


Pipe is a ubiquitous product in the industrial, commercial, and residential industries. It is 
fabricated from a wide variety of materials - steel, copper, cast iron, concrete, and various plastics 
such as ABS, PVC, CPVC, polyethylene, and polybutylene, among others. 


Pipes are identified by nominal or trade names that are proximately related to the actual 
diametral dimensions. It is common to identify pipes by inches using NPS or Nominal Pipe Size. 
Fortunately pipe size designation has been standardized. It is fabricated to nominal size with the 
outside diameter of a given size remaining constant while changing wall thickness is reflected in 
varying inside diameter. The outside diameter of sizes up to 12 inch NPS are fractionally larger than 
the stated nominal size. The outside diameter of sizes 14 inch NPS and larger are equal to the stated 
nominal size. Pipe wall thicknesses are specified by schedule number with 5 being the thinnest and 
160 be the thickest. An older designation scheme for pipe wall thickness which still enjoys popular 
usage indicates nominal weight. This labeling system is depicted in Table 1 on page 2. A more 
detailed dimensional table for NPS pipe is provided in Table 13 on page 77. 


When was It Used? 


Clay pipes have been found in excavations dated as early as 4000 BCE. They were used in 
Mesopotamia, the Indus Valley civilization, the Minoan civilization, and of course the Roman 
Empire (which also used lead pipes).' People have used pipes for thousands of years. Perhaps the 
first use was by ancient agriculturalists who diverted water from streams and rivers into their fields. 
Archeological evidence suggests that the Chinese used reed pipe for transporting water to desired 
locations as early as 2000 BCE.” 
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TABLE 1 Cross Reference of Pipe Wall Thickness Designations 
Schedule No. Weight Abbreviation Description 


5 = = 
10 LG Light gauge 
= LW Light weight 
20 = — 
30 = — 
40 ST Standard weight 
60 Em — 
80 XS Extra strong 
100 = — 
120 — — 
140 — — 
160 XX Double extra strong 


Early American settlers knew nothing of lead or iron pipe - they 
knew only to build with wood, the country's bounty. Water pipes 
were made of bored-out logs, felled from hemlock or elm trees. 


The trees were cut into seven-to-nine-foot lengths, with trunks 
being 9-10 inches in diameter. The interior was drilled or bored 
with steel augers. One end was rammed to form a conical shape, and logs were jammed together in 
series, using a bituminous-like pitch or tar to seal the joints. Sometimes the logs would be split and 
hollowed out, put together, and connected with iron hoops. A gravity water system would be set-up, 
starting from a spring or stream on high ground, allowing water to flow downhill to the house or 
farm. In the early 1700s, New York, as well as Boston, had constructed a wooden pipe system under 
the roads, and sold water at street pumps or hydrants. Wooden pipes were common until the early 
1800s when the increased pressure required to pump water into rapidly expanding city streets began 
to split the pipes. A change was made to iron. In 1804, Philadelphia earned the distinction as the 
first city in the world to adopt cast iron pipe for its water mains. It was also the first city in America 
to build large scale waterworks as it drew upon the ample supply of the Schuykill River.? 


Early pipe size selection was simple. The original pipe fabricators/layers of old were for the 
most part, not concerned with determining optimum pipe sizes since mere timber availability 
dictated the diameters installed. 
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Why is It Round? 


Aside from the use of the round tree by our predecessors due 
to its natural availability, pipe with a circular cross section is 
preferred over other geometric cross sections for a host of 
reasons. A quick mathematical examination reveals that a 
circular cross section results in the least surface area per unit 


of volume. Surface area is important because it relates 
directly to the amount of material required to fabricate the conduit, as well as the amount of 
protective coating or insulation required to cover it. The circular cross section is hydrodynamically 
more efficient than non-circular sections in that it presents less contact surface area for fluid friction 
for a given volumetric flow rate. Additionally, a curvilinear profile promotes a smooth flow path 
devoid of abrupt inside corners or pockets that promote localized fluid 
stagnation. Structurally, the curved shell produced by the circular pipe cross 
memerane A section allows a multiplicity of alternative stress paths and gives the optimum 


STRESS 


form for transmission of many different load types. This holds true for both 
internal loads like pressure and external loads such as those produced by back 


fill material or internal vacuum conditions. The single curvature allows for a 
simple fabrication process and is very efficient in resisting loads. Pressures are resisted very well by 
the in-plane behavior of a shell. Because this pressure difference is across a curved surface, a hoop 
of circumferential membrane stress is developed in the pipe wall. 


The Optimum Pipe Size 


Optimum pipe size denotes the best pipe size. From a simplistic 
standpoint, the best pipe size is obviously the smallest size that 
will accommodate the application at hand. From a realistic 
standpoint, optimum pipe size can have many meanings with 
proper consideration of the application. Optimum can mean 
economically efficient over the life of a system. Optimum can 
mean, as in the case of sulfuric acid, the pipe size which would 
limit the fluid velocity to a value which prevents pipe wall 


erosion in elbows which ultimately results in structural failure. 
Optimum can mean, as in the case of fluids with suspended solids, that pipe size which would 
produce a predetermined fluid velocity which is known to sustain the solids in suspension. 
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Too often, optimum pipe size is confused to be limited to mean most economic pipe size. 
Moreover, in addition to meaning satisfactory and maybe the most economical, optimum pipe size 
means that diameter which acts or produces the required effect with a minimum of waste and 
expense. This size can be far from the most economic due to specific process restraints. 


Quick Reference Chart of Useful Pipe Size Equations 


To gain immediate benefit from this course, this section provides a summary chart of formulas for 
determining approximate pipe sizes. The chart is presented at this point, before embarking on a 
detailed treatment of incompressible versus compressible fluid flow, single versus two phase flow 
patterns and the like. It offers quick access to formulas to estimate probable or target pipe sizes. 
This listing contains both rational and empirical equations that have evolved through industry's 
experience. They do not necessarily take into consideration possible mitigating parameters such as 
erosion, solids suspensions, or slug flow ramifications. These and other important considerations are 
covered later in the course. The mathematical expressions which follow are combinations of both 
simplified distillations of rational parent equations as well as rules-of-thumb. A list of the 
nomenclature used in the formulas and equations is provided immediately following the quick 
reference table. 


In the interest of brevity, equation nomenclature and units 
of measure are not provided separately with each equation. 
The student is therefore encouraged to frequently review (or 
better still. print) the nomenclature listine of Table 3. 


Units of measure for a specific equation variable are generally 
constant between equations. The exception is the flow variable “Q” 
which can be either cubic feet per second (ft'/sec) or gallons per 
minute (gal/min; gpm) depending on the need for unit consistency 
within the equation. Students are encouraged to conduct mental 
dimensional analysis when viewing equations and worked examples. 
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TABLE 2 Quick Reference Summary Chart & Collection of Pipe Size Formulae 
(Important: refer to Table 3 Equation Nomenclature and Units) 


EXPLANATION EQUATION RANGE & LIMITATIONS 


Pressurized Flow of Liquids 


Q in gal/min. 


0.2 
Darcy-Weisbach PEN | 0.0311 f LQ? | 
p h 
T 


Liquid general flow H Q : : 
EAR d= TA Q in gal/min. 
"M 1/3 
Nominal pipe size for 7 Q | : 
non-viscous flow gc 12) c 2 QUADR 
Nominal pipe size for E ; 
MR am E E = 02540 Q x 100 gal/min 
Pump suction size to 
limit frictional head d = 40.0744Q Q in gal/min 
loss 


frictional head loss 


Pressurized Flow of Gases 
Gas general flow d- 029 QT Q in standard ft’/min; T in °R; P in 
equation | 0 YA PV Ib/in? absolute. 

: 104 


Gas general flow Isothermal fully turbulent flow; Z in 
equation " miles; Q in standard ft*/day. 


Minimum diameter to 
limit erosional gas 
flow 


Weymouth gas flow » WGL Isothermal fully turbulent flow; Z in 
equation E 7868( gu P) miles. 


Two-Phase Flow 


Relief valve discharge _ | WNZ T/M, Isothermal flow of an ideal gas for 
flashing flow ~ | 408245MP selected Mach Number M. 
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Suitable pipe size for With known two-phase mixture 
gas/liquid flow velocity V. 


Steam & vapor 
general flow equation 


U Ww | 
| Manning formula for 7 0.375 eS 
uns d | Based on y/d = 0.983 ; Q in ft/sec. 

maximum flow 5 
Manning formula for Spee 
= . 1 3 
full pipe flow d z| Based on y = d ; Q in ft/sec. 


Minimum opening 
size of self-venting Qin gal/min; 2 Ed € 18; y/d — 0.75 
side entry overflow 
Minimum pipe size for 
self-venting vertical 
pipe flow 


Where r is ratio of annual flow area to 
pipe cross sectional area; Q in gal/min. 


TABLE 3 Equation Nomenclature and Units 


U.S. customary 


Symbol Definition units 

A Cross sectional flow area ft 

A Exponential correlation constant Dimensionless 

b Settling velocity constant Dimensionless 

C Total cost U.S. $ 

C Total life cycle cost present value U.S. $ 

C; Initial or first cost U.S. $ 

C. Annual operating cost U.S. $ 

E Mean volume fraction Decimal percent 

D Inside pipe diameter ft 

Dp Solid particle diameter in 

d Inside pipe diameter in 

FU Fixture Unit gal/min 

f Colebrook friction factor Dimensionless 
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TABLE 3 Equation Nomenclature and Units (continued) 


Symbol 


G 


WWE MSS ws eS Oe ee gn 


oN 


Vsu 


= SAST 


Definition 


Specific gravity 

Gravitational constant 

Energy input to overcome friction 
Annual cost escalation rate 
Froude number 

Critical static head 

Frictional head loss 

Settling velocity constant 
Length 

Mach number 

Molecular weight 

Hindered settling exponent 
Breakeven time 

Reynolds number 

Useful life of project alternative 
Manning friction factor 

Settling velocity constant 
Pressure 

Wetted perimeter 

Gas-phase pressure 
Liquid-phase pressure 
Volumetric flow rate 
Volumetric flow rate 

Hydraulic radius 

Ratio of annular flow areas 
Slope (hydraulic gradient) 
Saturated steam temperature 
Gas temperature 

System annual operation time 
Mean fluid velocity 

Mean fluid velocity 

Single particle settling velocity 
Hindered settling velocity 
Specific volume 

Saturated steam specific volume 
Weight flow rate 

Gas two phase weight flow rate 
Liquid two phase weight flow rate 
Steam quality 
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U.S. customary 
units 


Dimensionless 
ft/sec? 
Horsepower 
Decimal percent 
Dimensionless 
in 
ft of fluid 
Dimensionless 
fi 
Dimensionless 
Ibm /Ib-mol 
Dimensionless 
Years 
Dimensionless 
Years 
sec—ft 
Dimensionless 
Ib;/in? 
ft 
Ib;/in? 
Ib;/in? 
ft?/sec 
gal/min 
ft 
Dimensionless 
Dimensionless 
°F 
YR 
hr 
ft/sec 
ft/min 
ft/sec 
ft/sec 
ft?/Ib,., 
ft/Ib,, 
Ib, /hr 
Ib, /hr 
Ib, /hr 
Decimal fraction 


-1/3 
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TABLE 3 Equation Nomenclature and Units (continued) 


U.S. customary 


Symbol Definition units 
is Gas compressibility factor Dimensionless 
y Liquid flow depth ft 
Greek symbols 

a Piping system alternative choice Dimensionless 
p Piping system alternative choice Dimensionless 
^ Differential change Dimensionless 
E Pipe absolute roughness ft 
Nm Electric motor efficiency Decimal percent 
Np Pump hydraulic efficiency Decimal percent 
u Absolute viscosity cP 
v Kinematic viscosity ft/sec 

Fluid density Ibm /ft 
pr Mixed phase fluid density Ib,, /fe 
pr Solid particle density Ib,, /ft? 
X Mole fraction Decimal fraction 


When units other than those listed are used, it will be expressly noted 
A at the time of equation presentation. 


Pipe Flow Basics 
Flow Regimes 


During the 1880s an Irish Engineer and physicist named Osbourne Reynolds conducted experiments 
by visually monitoring flow patterns in glass tubes with dye injected fluids. His observations 
resulted in the now famous dimensionless quantity which bears his name, the Reynolds number, 
Nre. The Reynolds number relates the physical and geometric properties of fluid density, velocity, 
viscosity, and pipe diameter. His studies showed that in essence three flow regimes exist: 


Laminar (smooth) for Nre < 2,000; 
Transitional for 2,000 < Nre < 4,000; and 
Turbulent for Nre > 4,000 
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When colored filaments are injected into a flowing fluid stream, their appearance assumes 
the configurations shown in the diagram below for each of the three regimes. 


A LAMINAR FLOW 
= 

A TRANSITION FLOW 
=h (CRITICAL ZONE) 
LILL ey 

ey TURBULENT FLOW 


Figure 1 — Flow Regimes 


Flow Variables 


Part of the difficulty in working with fluid flow problems is the wide array of different variables 
that come into play. It is helpful to categorize these variables. The fluid dependent variables are: 


e Fluid density; 
e Fluid viscosity. 


The system dependent variables are: 


e Surface roughness of pipe interior; 
e Length of flow path. 


The basic flow parameters are: 


e Flow velocity; 

e Pipe size; 

e Flow rate; 

e Pressure loss due to frictional drag. 
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The fluid dependent variables remain essentially constant for most problems of interest. 
Similarly, the roughness is constant. Therefore, this leaves pressure drop and the four variables that 
have been characterized as the basic flow parameters, that must be balanced to select pipe size 
properly.^ 

Invariably, there are four basic equations which are used in the solution of pressurized fluid 
flow problems: 

Continuity equation modified for volumetric flow; 
Reynolds number; 

Colebrook friction factor and; 

the Darcy-Weisbach pressure drop. 


Each of these will be used throughout the course. If necessary, the student should review 
PDHcenter.com Course No. M212 for a detailed treatment of the first two items above. 


Fluid Friction 


It is impossible to study pipe size selection without an understanding of fluid friction and in 
particular the fluid friction factor f. A great deal of technical literature is available, in print and on 
the Internet, on the subject of computing system friction. Any of this information can be 
successfully employed depending upon the student's preference, past experience, and applicability 
to the system and liquid being considered. In all cases, the line losses vary directly as a function of 
the square of the mean fluid velocity. 

The friction factor 1s a critical variable regardless of the technique employed to determine 
frictional head loss. It is therefore important to have a background knowledge of its history and 
development. In 1938 British Engineer C.F. Colebrook demonstrated that for Nre > 3,000, a fluid 
friction factor existed which was a function of both the value of Nre and the relative roughness (¢/D) 
of the pipe. No approximations here; the Colebrook factor is extremely accurate. The Colebrook 
friction factor is defined as 


1 g€ 251 
| [1] 


A = [sss e 


This hideous looking mathematical creature is implicit in f which requires an iterative 
process for convergence to a solution. Going forward, a host of later experimental Engineers 
developed approximations of this friction factor concept that were explicit in f, culminating in a 
graphical representation of a friction factor relationship produced by L.F. Moody in 1944. A 
schematic diagram of what has become to be known as the Moody chart is shown in Figure 2. It is a 
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collection of plots of friction factors versus Reynolds number for a variety of relative roughness 


values (&/D). 


0.052 
D. 04 — 


D. m | 
0.036 


0.032 


0.026 


f 


0.024 


0.020 


Friction Factor 


0.016 


0.012 


LIII " 
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ES SEE EHH 
SL 
N —— an 
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oa Msn 
[ — | [ 1 [I| 
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f 4 = = qmm = 


10 2 345 w 2 un ui 10 2 345 10" 
Reynolds aan 


Figure 2 — Moody Chart or Diagram 


The Moody chart was used extensively prior to the advent of the programmable calculator 


and the desktop computer. Present day desktop computing power has eliminated the tedious 


exercise of manually solving the Colebrook equation by iterative convergence necessitated by the 


fact that the Colebrook equation is implicit in f. 


The drudgery (and resulting approximation) associated 
with the use of the Moody chart to determine friction 


factors has been gratefully replaced with any of a number 
of free applets available on the Internet. 
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Alternative methods and approximations for determining the value of f have been studied 
ad infinitum and the number of empirical relationships associated with it abound. No attempt will 
made in this course to list or describe these various formulae or recommend one particular formula 
over another. For no particular reason, when the need for a value of f is required, the Colebrook 
value will be referenced. 

While the Moody chart provides a quick method of f determination, an extension of a 
concept derived from the Moody chart can be used to rid us of the burden of evaluating fat each 
flow rate. The Moody chart f values become asymptotic (f~ constant) in the zone of complete 
turbulence. With a close examination of Figure 2, and using water as the fluid, a fully turbulent 
fluid flow friction factor can be shown to approximately vary with pipe diameter in accordance with 
Table 4. 


TABLE 4 Approximate Asymptotic Moody Friction Factors for 
Clean Commercial Steel Pipe with Flow in Zone of 
Complete Turbulence 


Nominal Pipe Size (inches) Friction Factor (f) 
1 0.023 
2 0.019 
3 0.018 
4 0.017 
6 0.015 
8-10 0.014 
12-16 0.013 
18-24 0.012 
30-48 0.011 


The Fanning friction factor can be cause for much 
consternation. Some Nz. versus &D fluid friction charts list 
values of the Fanning friction factor along the vertical axis in 


lieu of the Colebrook friction factor. It is important to 
remember that the Fanning factor = Colebrook + 4. 
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Frictional Head Loss 


If a universal equation exists with regard to the flow of fluids it would have to be the Darcy- 
Weisbach formula. Within limits, in can be used to analyze both incompressible and compressible 
fluids. All of its variables, with the exception of f, can be rationally derived. 


peg [2] 


Without derivation, the Darcy-Weisbach energy loss equation can be rewritten in terms of 
pipe diameter d with Q in gallons per minute, to appear as 


2 0.2 
2 p| Sree T 
f 


This form allows for the direct determination of an optimum pipe size. 


A Trade-Off 


One can always trade-off pipe size against velocity and pressure drop. Greater flow for a given pipe 
size can be obtained if a higher pressure drop is acceptable. Given the fact that static heads are 
constant regardless of pipe size, comparable pressurized pumping costs are purely attributable to 
fluid friction. This leads us to the subject of selecting the most economic pipe size. 


Economics 


Comparative Analysis 


Simply stated, there is no easy approach to determine an economic pipe size for pressurized 
systems. There have been volumes written on optimum economic pipe size. Any method which only 
looks at flow rate to determine pipe size (or vice versa) cannot cover all the cases. From a purely 
economic standpoint involving pressurized flow, the optimum pipe size is the diameter that 
minimizes the life cycle cost of a piping system. The first (installation) cost and the recurring costs 
associated with operation and maintenance are the components that are examined using the concept 
of time-value of money. It follows that there is a compromise between large diameter pipe (higher 
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first cost, lower recurring costs) and small diameter pipe (lower first cost, higher recurring costs) for 
a given service. The only tried and true method involves detailed comparative analysis. 


There is no shortage of coverage in the open technical 
literature on the subject of picking directly the most 
economic pipe size. Practically all of the methods presented 
depend on date sensitive cost data which quickly becomes 


archaic after publication. Simplified methods of 
determining economic pipe size such as nomographs which 
incorporate material, energy, and labor costs, are 
problematic for the same reason. 


In order to illustrate fully the concept of hardware versus pumping cost, an example will be 
provided which introduces the comparative cost method known as break-even analysis. 


Break-Even Analysis 


Many capital investment alternatives cannot be evaluated by conventional means of worth. Such a 
case is one in which neither of the presented alternatives can be justified solely on the basis of an 
inherently derived savings. Break-even analysis becomes useful in such situations. For the same 
pressurized fluid flow application there is a tradeoff between a large pipe size (higher first cost C; 
with lower operating cost C,) and small pipe size (lower first cost C; with higher operating cost C,); 
consequently, from a total cost standpoint, there exists a point during the life of the project where 
the attributes of each are equal. This is known as the break-even point. When inherent savings are 
not applicable, it can be useful to identify the periods that each alternative is economic-ally 
attractive during the course of the installation's expected life, n. Piping systems often have a life 
spans of 20+ years. 


The concept of break-even analysis uses the following thought process: 
The total cost of one alternative, let's call it a, can be expressed as 
Ca = Cai to, [4] 


The total cost of the second alternative f could be indicated by 
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Cy = Cg; + n Cg, [5] 


To express the relationship in terms of break-even point of these alternatives, let the total cost of 
each be equal: 
CutsG,, = C5; + Cg, [6] 


Rearranging this equation in terms of n results in a relationship in which the break-even point (in 
time) is indicated: 


n 


e = L7] 
C Cao 


a 


This relationship is shown graphically in Figure 3. 


«s 
advantage 


Cost ($) 


B’s 


1st cost advantage 


Figure 3 — Break-Even Point of Two Alternative Systems 


e Time 


Studying this diagram shows that the break-even point occurs where the total costs for the 
individual alternatives are equal. Put another way, the break-even point is where the total cost of 
piping system f (which represents the lowest first cost but has the higher annual operating expense) 
has reached the value of the total cost of piping system a (higher first cost but lower annual 
operating cost). Let us look at an example of break-even analysis for a hypothetical pressurized pipe 
system which offers two alternatives. 
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It should be noted that annual operation time can be very 
critical in the determination of optimum pipe sizes for 
externally pressurized flow. It is intuitively obvious that 
operating cost is of little significance where the annual duty 


cycle (pumping or compressing hours) is very low. In these 
systems the smallest pipe size otherwise allowable is the most 
economical and therefore the optimum. 


Example 1 Evaluation of Two Alternative Piping Systems 


Problem: A choice exists between investing in piping system a or piping system f. Alternative a 
utilizes large diameter piping with an initial investment of $100,000 and an annual operating energy 
expense of $20,000. Alternative f employs smaller diameter piping resulting in an initial outlay of 
only $80,000 but due to the increased pumping cost, is estimated to have an annual operating 
expense of $25,000. Both systems have been estimated to have a useful life of 15 years. From an 
economic standpoint, which alternative system should be installed? 


Given: Cai = $100,000 ; Cao = $20,000 ; Cj; = $80,000 ; C5, = $25,000 ; n = 15 years. 


Find: The break-even point and the most economical alternative system. 


Solution: 
Using Equation [7] and substituting the known values results in a break-even point of 


Cu - Cpi _ 100,000- 80,000 _ , 
= = = ar 
Am 25000-20000 ~ 7 


n 


ao 


Conclusion: Since the useful life of both systems is indicated at 15 years, piping system a is 
decidedly more attractive because of its lower yearly operating energy expense. Were the useful life 
of the two systems to be less than 4 years, alternative / would be the most attractive. 


The present value of the total cost of each alternative in Example 1, using the time value of 
money, would need to be determined in order to properly evaluate the comparative costs. This leads 
to the concept of Life Cycle Cost Analysis. 
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Life Cycle Cost Analysis 


The life cycle cost (LCC) of any piping system is the total /ifetime cost to purchase, install, operate, 
and maintain that piping system. LCC analysis requires the evaluation of alternative systems based 
on the present value of these monetary outlays. 

Piping system project scale must be considered when allocating time spent in the life cycle 
cost analysis exercise. It stands to reason that time expended should be proportional to the physical 
size of the piping system. To wit, the design of an overland gas transmission pipe line would be 
given heightened life cycle cost analysis emphasis compared to a simple in-plant piping system. A 
piping system of relatively short length may not warrant any consideration whatsoever in this 
respect. 

Equivalent present value analysis requires the conversion of all future cash outflows 
(recurring annual costs) to the present. As such, it requires the consideration of time value of 
money using a term known as the uniform series present value factor. Comparison of the equivalent 
value of competing alternatives allows the choice of the most desirable alternative on the basis of 
economics. 

When the estimated useful project life is identical for each alternative, the alternative with 
the lowest present net value should be selected. Recall from engineering economics that the uniform 
series present value factor is, 

(14 i)'-1 
dae» B 
The annual operation cost variable C, must be modified by this factor in order to arrive at 


the total life cycle cost present value: 


(14 i)*-1 
C= GtC, iti)” [9] 
Here is a step-wise approach for determining the most economic pipe size: 

1. Knowing the required process flow rate, select a trial pipe size from an appropriate equation 
listed in Table 2; 

2. Using an accurate piping arrangement drawing, as opposed to a piping and instrument 
diagram (P&ID), determine the exact physical length, including equivalent lengths of the 
various valves, fittings, and equipment, of the system. Put another way, identify and 
quantify all components of the system which result in flow resistance; 
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3. With the given fluid properties, required flow rate, and pipe material and length, calculate 
the Reynolds number, relative roughness, friction factor, and head loss due to friction; 

4. Assume appropriate equipment operating efficiencies. For pumps use maybe 6096; for 
electric motors use 90%. Determine the required energy input to overcome system fluid 
friction losses; 

5. Basedon the stated running hours, compute the system's annual operating energy cost based 
on the energy input calculated in Step 4; 

6. Obtain a manufacturer's price quotation for the pump and driver based on the required flow 
rate and the calculated head from Step 3; 

Assume that annual maintenance costs will be 4% of the price obtained in Step 6; 

8. Sum the annual costs of energy and maintenance to arrive at the variable C,, total annual 
recurring operating cost. (Include any other recurring costs that might be applicable here); 

9. From the detailed piping arrangement, estimate the material and labor cost to install the trial 
pipe size system. Sum this cost with the equipment price quotation of Step 6 to arrive at the 
variable C, the initial capital cost; 

10. Use Equation [9] to determine the total life cycle cost present value for the trial pipe size; 

11. Repeat Steps 2 through 10 for one nominal pipe size below and one nominal pipe size above 
the trial pipe size selected in Step 1; 

12. Compare the three separate alternative present values and select the pipe size corresponding 
to lowest cost. 


An illustration of this step-wise economic analysis is shown in the worked example which 
follows. This example is abbreviated and only the details of the initial iteration cycle are presented; 
however, the results of the entire analysis are summarized in Table 5. 


Example 2 Selection of Most Economic Pipe Size for a Given Process Condition 


Problem: ^ clean fluid is to be continuously pumped at a rate of 200 gallons per minute through 
960 feet of schedule 40 alloy steel pipe. The fluid has a specific gravity of 0.8 and a viscosity of 1.1 
cP at the pumping temperature. There is no flow independent head, i.e. static head, in the process. If 
annual cost escalation is taken as 5%, annual maintenance cost is 4% of capital equipment cost, and 
energy cost is 7¢ per kwh, what is the most economic pipe size for this operation? The system is 
expected to have a useful life of 15 years. 


Given: t,— 8,760 hours ; Q — 200 gpm ; L= 960 feet of steel pipe (e = 0.00015 feet) ; G=0.8 ; 
4L — 1.1 cP ; i=0.05 ; maintenance allowance = 0.04 ; power cost is $0.07 per kwh ; n= 15 yr. 
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Find: The most economic pipe size d based on lowest life cycle cost present value of the system. 


Methodology: The step-wise cost analysis approach just described will be used. In the interest of 
brevity, for this example many of the specific calculation details will not be presented. 


Solution: 
O For Q = 200 gpm, from the nominal pipe size equation in Table 2 determine that d = 3.503 


inches. From Table 13 on Page 77 select a 3 inch NPS schedule 40 pipe with d — 3.068 inches as a 
trial size. Later, for comparison, examine 2% inch NPS and 4 inch NPS as economic alternatives. 


@ With the given fluid properties, required flow rate, and pipe material and length, calculate, 


V = 8.7 ft/sec 

Nre = 149,694 

&D - 0.00059 

f = 0.01976 

hy = 87 feet of fluid 


Q Recall that input power to a pump is given by the formula presented below. Using the assumed 
pump and motor operating efficiencies offered in the course, determine the required horsepower to 
overcome the system friction in the 3 inch NPS schedule 40 pipe to be, 


- Qh; G (2008708) —— 651 hp 
39605,7, (3960)0.60)(0.90) 


® The annual operating cost for this system is, 


_ ($0.07) kwh  6.5lhp 8760hr 
^ — kwh 134lhp-hr 1 1 


- $2,977 [round-up to $3,000] 


© Obtain a manufacturer's quotation for the pump/motor set. For this iteration (probably 72 
horsepower and 90 feet of head) assume this value to be $13,000. From this equipment cost we can 
estimate an annual maintenance cost to be (0.04)($13,000) = $520. The value of the recurring 
annual cost C, in Step & becomes C, = $3,000 + $520 = $3,520. 


© From a detailed piping arrangement of the proposed system, estimate the material and labor cost 
to install the trial 3 inch NPS pipe system. For this example, assume this cost is $21,000. The initial 
capital cost is the sum of the piping and equipment cost, or C; = $13,000 + $21,000 = $34,000. 
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© Total life cycle cost for the 3 inch NPS size system consists of the sum of the first cost and the 
present net value of the recurring annual costs: 


(i+i)"-1 (1+ 0.05)'°- 1 


C=C ae = $34,000 + $3,520. — 15 
i(1+ 0.05) 


i es mp - $70,540 
i(1l4 i) 


When Steps © through @ are repeated for 2% inch and 4 inch NPS schedule 40 pipe, the 
summarized data in Table 5 is generated. 


TABLE 5 Summary of Economic Variables for Three 
Alternative Pipe Systems of Example 2 


Nominal Pipe Size (NPS) 

Parameter (variable, units) 2% inch 3 inch 4 inch 
Diameter ( d, inches) 2.469 3.068 4.026 
Mean fluid velocity (V, ft/sec) 13.4 8.7 5.0 

Reynolds number (Nr) 186,012 149,694 114,074 
Relative roughness (¢/D ) 0.00073 0.00059 0.00045 
Colebrook friction factor (f) 0.02001 0.01976 0.01975 
Head loss due to friction (A, feet) 260 87 22 
Power input (HP, hp) 19.45 6.51 1.65 

Annual recurring costs (C,, $) $ 9,120 $ 3,20 $ 1,150 
Initial capital cost (C;, $) $ 30,000 $34,000 $ 38,000 


Life cycle cost present value (C, $) $124,660 $70,540 $ 49,940 


Q The alternative with the lowest present net value is the 4 inch NPS. Therefore it is the most 
economic pipe size. 


The content just presented was limited to pressurized piping systems. Gravity flow does not 
involve external energy input and therefore LCC analysis with regard to power consumption is 
inconsequential. Gravity flow design applications are often quite challenging however. Coverage of 
the topic of gravity flow also provides the opportunity to introduce the concept of flow patterns 
which, by-the-way, are manifested in both gravity and pressurized piping systems. 
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Gravity Flow 


Introduction 


Some of the trickiest situations concerning pipe line sizing involve flow solely under the influence 
of gravity. This is true in both single-phase (full pipe) and two-phase (partially filled) flow patterns. 
Two-phase flow occurs in a pipe when a gas and liquid are coexistent. This phenomenon is 
characterized by flow patterns. 


Flow Patterns 


Flow pattern is distinguished from flow regime where the fluid is said to be laminar or turbulent. 
Depending on the gas to liquid ratio and the mixture velocity, flow patterns range from bubble flow 


to dispersed flow. See Figure 4. 


Ae O ee 
Ae O DON diii 


dir e dui d 
= STRATIFIED-WAVY FLOW 


DISPERSED FLOW 


(LIQUID ENTRAINED 4S 
SPRAY IN GAS PHASE) 


Figure 4 — Flow Patterns in Horizontal Pipes 
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Flow patterns can change with pipe orientation alone. See Figure 5. 


i 


VERTICAL BUBBLE VERTICAL SLUG 


FLOW 


‘ |] 
FLOW ANNULAR FLOW 


Figure 5 - Flow Patterns in Vertical Pipes 


Plumbing Considerations 


DWY PLUMBING STACK 


Selecting the Optimum Pipe Size 


Obviously the fluid mechanics are the same in process piping and 
plumbing; however, the terminology and the design approach to 
select pipe size can be quite different. For instance, vertical runs of 
processing piping have no specific name; vertical runs of plumbing 
waste or drainage lines are often called stacks, while supply lines 
are known as risers. The constituents in the gas phase of two phase 
process piping flow are varied. The countercurrent core flow of gas 
in plumbing annular flow is essentially air. Venting in plumbing 
stacks is needed to preserve the trap, that liquid seal which 
prevents the introduction of obnoxious odiferous vapors into 
occupied spaces. Gravity flow process piping is not necessarily like 
plumbing applications where the optimum pipe size must be 
selected to produce a low velocity liquid flowing vertically 
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downward as an annular film with displaced gas being allowed to flow counter-currently (vent 
stack). (The term DWV means drain/waste/vent). This is not true in industrial process piping where 
the gas stream contained within the core is often times flowing concurrently, not counter- 
concurrently, with the liquid annular film. 

Building designers and building mechanical engineers use a simplified design procedure 
from building codes to select drainage pipe size within buildings. It is customary to size plumbing 
using a cookbook approach. Building codes generally dictate minimum pipe sizes based on the 
fixture unit concept. The procedure assigns fixture units (FUs) to various waste generating 
plumbing fixtures, such as sinks and showers, then establishes authoritatively the permissible pipe 
sizes and slopes for a given number of total fixture units contributing to a point in a drainage 
system. The designer determines the total number of FUs and selects an accommodating pipe size. 
Table 6 is a typical example of just such a code listing. 


TABLE 6 Typical Plumbing Code Fixture Unit Loading Chart 


Pipe Size (inches) Maximum Fixture Units 
1% 1 
1% 1 
2 8 
2⁄2 14 
3 35 
4 216 
6 720 
8 2640 
10 4680 
12 8200 


Stack flow has been studied extensively? because of the deleterious effect which can be 
created when critical flow rates are exceeded. It has been found that slugs of water and the resultant 
violent pressure fluctuations do not occur until the stack flow exceeds one-third full. Plumbing 
stacks are generally sized to insure that flow corresponds to r < '^ full because higher flows result 
in pressure fluctuations which compromise the trap seal. The maximum permissible flow rates in 
the stack to prevent slug flow can be expressed by the formula: 
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4A 
Q = 278r? qus where r= EE) [10] 


The variable r should not be confused with 4A d 


hydraulic radius. It is the ratio of the cross sectional area r= ad? 
of the annular flow to the total cross sectional area of the 
pipe. See Figure 6. Equation [10] rewritten in terms of d 


to allow for a direct solution of pipe size yields, 


Q 3/8 
4 — un 


Il 
AS 


Figure 6 — Plan View of Stack Flow 


With regard to gravity flow in partially filled pipes, the 
term full is used somewhat inconsistently. In vertical pipes 
and stacks, 1⁄4 full and % full generally refers to the value of 
r, the ratio of flow area to inside pipe area. In horizontal 
and near-horizontal lines, % full and 7^ full generally refer 
to the value of y, the relative flow depth (see below). 


Figure 7 — Relative Flow Depth in Horizontal and Near-Horizontal Pipes 


Selecting the Optimum Pipe Size 24 
O 2008, 2015 Randall W. Whitesides, CPE, PE 


www.PDHcenter.com PDH Course M270 www.PDHonline.org 


After about 30 feet, gravity flow liquid reaches terminal velocity in vertical runs of pipe. 
The flow rate in Equation [11] is in gallons per minute once terminal velocity is achieved. The 
formulas developed for terminal velocity and terminal length, without derivation? 5, are given in 
Figure 8. 


Terminal Fluid Velocity 
v= 3(oa)* 


reached at 7 = 0.052 y? 
from source 


HYDRAULIC JUMP 


Figure 8 — Vertical to Horizontal Transition Flow 


At the transition of vertical to horizontal flow at the base of the pipe, the fluid velocity 
decreases after turning horizontally. This change in kinetic energy results in a hydraulic jump. 
While the hydraulic jump is normally of little consequence, the fluid characteristics upstream of the 
jump can essentially close-off the vapor phase flow, ultimately resulting in potentially damaging 
slug/plug flow pulsations. 
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Example 3 Selection of Minimum Pipe Size for a Given Waste Flow Condition 


Problem: A proposed high-rise building project's design specifications indicate that plumbing 
stacks shall be sized to flow no more than "4 full to prevent damaging slug flow. What diameter 
pipe should be used to permit waste water to flow at the rate of 550 gallons per minute? Determine 
the terminal fluid velocity, the vertical distance from the 550 gpm flow source at which this velocity 
will be achieved, and the approximate lateral distance from the stack at which a hydraulic jump will 
form after the DWV line turns to a horizontal orientation. 


Given: r= '4 (0.25) [assumed, see below] ; Q = 550 gpm. 


Find: A target pipe diameter d, terminal velocity V, terminal length L, and a horizontal distance to 
jump formation. 


Solution: 


® While it is not stated in the problem, as is customary practice, 1⁄4 full shall be assumed to mean 
the ratio of flow area to total pipe area. Using Equation [11] and the given quantities, an inside pipe 
diameter is calculated to be, 


3/8 3/8 
Q 550 


Q From Table 13 on page 77 select a 8 inch NPS schedule 40 pipe with d= 7.981 inches 
Q From Figure 8, calculate the terminal velocity, vertical distance, and lateral displacement of the 
resulting hydraulic jump to be, 


2 0.4 l 550 y 

= = = ——— = 16.3 ft/ 

r (£ 3 7981 - 
L = 0.052V? = (0.052)(163)? = 13.8 feet 


10d = (10)7.981)12 = 6% feet 


Selecting the Optimum Pipe Size 26 
© 2008, 2015 Randall W. Whitesides, CPE, PE 


www.PDHcenter.com PDH Course M270 www.PDHonline.org 


Civil Engineering Applications 


Open channel flows are modeled most frequently by the Manning formula which in some special 
cases can also be used to select pipe size. Irish Engineer Robert Manning is attributed with the 
development of this simple, dimensionally homogeneous formula for gravity flow. The Manning 
formula’ is semi-empirical. In its original form it is based in terms of fluid velocity, 


V= : R23 gi [12] 


where V is the liquid cross-section mean velocity and R is the hydraulic radius. Hydraulic radius 
should not be confused to mean half of the pipe inside diameter. Hydraulic radius is, 


A 


R= 
wetted perimeter us] 


The Manning formula does not account for entrance or exit losses, bend losses, losses due to 
flow through a partially open/closed valve, loss through reducers or any other minor losses. It is 
best used for uniform steady state flows. Uniform means that the pipe size remains constant and 
steady state means that the velocity, discharge, and depth do not change with time. These 
assumptions are rarely ever strictly achieved in process conditions. 


The formulas, equations, and inequalities presented for gravity 
flow are, for the most part, empirical for water at 68? F (u — 1 cP). 
Stated reductions in flow rate of 10% have been noted in Reference 


9 for fluids with viscosity on the order of 10 cP. Consequently, the 
information presented must be used with due consideration when 
very viscous fluids are encountered. 


Experience has shown that the maximum gravity water flow does not occur when a near- 
horizontal pipe is flowing full. Fluid friction reduces flow rate when the ratio of depth to diameter is 
greater than y /d = 0.938. A rational mathematical solution verifies that the additional flow area 
gained by increasing the ratio's value greater than 0.938 is differentially small compared to the 
large additional pipe surface area presented for fluid friction. In a similar fashion, maximum fluid 
velocity can be shown to occur at y/d = 0.81. Armed with this knowledge, the standard form of 
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the Manning formula presented in Equation [12] can be rewritten to allow the determination of a 
pipe size which offers maximum gravity flow for a given flow rate and energy gradient: 


On 0.375 
= = 14 
d 1525 e [14] 


Unlike civil engineering applications of open channel flow where it is often difficult to 
accurately estimate the roughness coefficient n of a natural channel bed, the Manning roughness 
coefficients for commercially available pipe materials is fairly well established, varying within a 
range. These are shown in Table 7. Another useful rewrite of the standard Manning formula is one 
in which the pipe is flowing full’: 


0.375 
d= 1639 e [15] 


TABLE 7 Ranges of Manning's n 


Pipe Material Manning's n 
Cement 0.011 — 0.017 
Clay 0.010 — 0.015 
Ductile Iron 0.011 — 0.016 
Plastic 0.010 — 0.015 
Steel 0.010 — 0.014 


Example 4 Optimum Pipe Size for Near-Horizontal Gravity Flow Transfer Pipe 
Problem: Cool beer wort (whose physical properties are very similar to H20) is to be transferred at 
the rate of 600 gallons per minute between two brewery vessels through a near-horizontal stainless 
steel transfer line. The liquid levels in the two vessels will be automatically maintained to establish 
uniform, steady flow in this crossover pipe. Spatial limitations dictate that the maximum slope 
(gradient) for this line be 1:14. Discounting vessel liquid level effects at the pipe entrance and exit, 
and without concern for gas entrainment, what is the smallest pipe inside diameter which will 
accomplish this clarified wort transfer? 
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Given: n — 0.013 (stainless steel) ; O = 600 gpm = (1.34 ft?/sec) ; S = 1:14 (or 0.07). 


Find: The optimum (smallest) pipe diameter d 
for the stated process conditions. 


Methodology: Since no mention of 
pressurization is made, the decanting operation 
described implies a gravity flow partially filled 
pipe scenario. Accordingly, the Manning 
formula appears to be acceptable. The smallest 
pipe for this application will be the one which 
provides the maximum partially filled pipe flow 
depth for uniform 600 gpm steady flow. 


Solution: 


O Using Equation [10] and the given quantities, calculate a trial diameter of 


or 0.375 | EI 0.375 
d = |1525 “= = | 1525 -———— = 5.63 h 
| Js J0.07 ilic 


@ Find the nearest standard pipe size corresponding to the trial size found in Step (D. From Table 
13 on page 77 selecta 6 inch NPS schedule 10 standard pipe size with d = 6.357 inches 


Process Drainage 

Basically, gravity liquid flow through pipe from relatively large source volumes can be approached 
from three distinct design perspectives. These are (1) full pipe single phase flow; (2) gas entrained 
flow; and (3) self-venting flow. Of the three conditions, the full single phase flow design generally 
results in the smallest pipe size but it can be the most difficult to maintain. Problems arise when 
actual flow conditions do not match the design approach used. When such pipes do not run full, 
considerable amounts of gas can be drawn down by the liquid. The amount of gas entrained is a 
function of pipe diameter, pipe length, and liquid flow rate. 


Entrainment inhibits liquid gravity flow by raising 
the pressure drop through the piping above that for 


single-phase flow, and reduces the static head 
available for over-coming the pressure drop. 
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determined on the basis of full pipe flow but in fact the situation is entraining gas because a critical 
minimum head fails to be maintained. In this case an absorption column is being drained through a 


vertical pipe whose size has been designed for full flow. Absorption columns are vessels which are 


used throughout the petrochemical industry to accomplish mass transfers of fluids and gases. 


TREATED 


GAS FLOW 


OUT 


Liquid flow in ex- 


trained in outlet. 
The vessel level 
rises. 


falls rapidly. 


Figure 9 — Cyclic Surging in Gas Entrained Drain (adapted from Reference 9) 
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The higher liquid 
level reduces 
entrainment; but 
drainage is still 
restricted. Level 
continues to 
rise. 


! 


(5) Drainage flow -4—— 
remains high. 
Built-up gas in 
inlet line is releas- 
ed in a surge. En- 
trainment begins 
again. Conditions 
are —» 


! 


G) Higher level 4—— 
stops entrain- 
ment. Bubbles trap- 
ped in outlet con- 

tinue to restrict 

flow. Level rises, 

gas flow begins Y 
to be restricted. 


The original 
conditions have 
been reestablished 
and the cycle re- 
peats. 1 
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Froude Number 


William Froude was a 19" century British Engineer, hydrodynamicist, and naval architect. He 
formulated reliable laws for the prediction of ship stability and for the resistance water offers to 
marine craft. He is credited with the development of the dimensionless number presented here. He 
developed a dimensionless number that compares inertial and gravitational forces. 

In the context of this course, the Froude number is a dimensionless liquid flow rate term 
which can be used as a superficial volumetric comparator: 


y 
d [16] 
129 


Where V is based on the assumption of full-pipe flow.'° 


The number is useful in sizing gravity flow piping because extensive observations have 
yielded empirical relationships which define the interdependency of Froude number, static heads, 
flow rates, and pipe size. 


Pipe Size for Flooded Gravity Flow 

For full flow, the pipe size must be based on single-phase criteria. To avoid gas entrainment in the 
full-pipe-flow design, the motive force, i.e. static head, must be great enough to maintain a flooded 
inlet condition. According to Hills?, two simultaneous conditions must be met in order to insure full 
pipe single phase flow free from gas entrainment. These conditions which involve static liquid 
heights and Froude number values are given by the following inequalities: 


For full pipe single phase liquid flow from bottom outlets: 


2 0.25 2 
h > 10704 (22 add FS 16(4) [17] 
& 


Where Q is in ft/sec. 
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For full pipe single phase liquid flow from side outlets: 


2 
h > 201803 A and Res [18] 
g 


The reference datums for these relationships are shown in Figure 10. 


Figure 10 — Critical Static Head for Full Pipe Single-Phase Liquid Drainage 
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Pipe Size for Self-Venting Flow 


Self-venting drainage is generally carried out by means of a side outlet arrangement. Hills? suggests 
that two simultaneous conditions must be met in order to insure side outlet self-vented flow. These 
conditions, which involve a diameter at the entrance boundary to the outlet and the Froude number, 
are given by the following inequalities: 


04 
4 
d > 12/42] with F, < 03 [19] 
03mg 
The above inequalities insure that the depth of flow in the pipe will be less than half full at the 


entrance to the outlet. To insure self-vented gravity flow, a critical far field (away from the outlet) 
static fluid height less than 0.8d must be maintained. This fluid height datum is shown in Figure 11. 


Figure 11 — Static Head to Promote Self-Venting Drainage 
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A regression analysis by Whitesides of Hills? pipe flow data resulted in the following flow 
capacity correlation equation which related nominal pipe size to side outlet gravity flow for the 
arrangement shown in Figure 11: 


Q-123d?  [fo2«d < 18] [20] 


Where Q is in gallons per minute 


Written in terms of d, Equation [20] becomes, 


d = 0920"  [for5« Q < 1,600] [21] 


For a desired flow rate, a diameter for the outlet's entrance can be estimated directly from 
the Inequality [19], from which a fluid velocity and Froude number can be established. If the 
Froude number is less than 0.3, then it can be assumed that the selected opening size is acceptable. 
The design procedure for self-venting gravity flow is shown in Example 5. 


Example 5 Optimum Pipe Size for Self-Venting Gravity Flow from Vessel 


Problem: It is desired to vent liquid from an atmospheric tank at a minimum rate of 120 gallons per 
minute through a side outlet self-venting drain (overflow arrangement Figure 11). What size free 
opening should be made in the tank wall? What far field (away from the outlet's entrance) 
maximum liquid level should be allowed above the outlet's invert to avoid air entrainment? 


Given: Q = 120 gpm = 0.267 ft/sec 


Find: An acceptable opening diameter d and corresponding critical static head A which should not 
be exceeded to insure self-venting drainage without surging. 


Methodology: Using the Inequality [19], estimate an initial opening size. Knowing the required 
flow rate, calculate a corresponding fluid velocity for this diameter. Test Inequality [19] to 
determine if the Froude number is acceptable for this size. Repeat this process, increasing the size, 
until an acceptable size is found. 
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Solution: 


O Determine a trial pipe size from Inequality [19] for a self-vented side outlet configuration, 


04 
4Q (4)(0.267) . 
d > HER > 2 Or um > 6.30 inches 


@ Find the nearest standard pipe size corresponding to the trial size found in Step (D. From Table 
13 on page 77 select a 6 inch NPS schedule 10 standard pipe size with d = 6.357 inches 


®© With the required flow rate, determine the fluid velocity for the selected standard pipe size 
assuming full pipe flow. Use a rearrangement of the continuity equation. 


4 4)(0.267 
Q = o = (4X ) = 1.21 ft/sec 
A amd 7(6.357 / 12) 


ys 


® Calculate the Froude number for the selected standard pipe size and the computed fluid velocity, 


y 121 
F. = = = 
| d 4(6.357 /12)(32.17) 
12 


F, = 0.29 «0.3 .. by Inequality [19], the selected 6 inch NPS Schedule 10 pipe is o.k. for side 
outlet self-vented flow 


@ Maximum static head to insure 1⁄2 flow depth in the outlet will be 4 = (0.8)(6.357) = 5 inches 
(& Use Equation [20] to determine the probable flow rate for the selected drain outlet pipe size, 


Q-1234?? = (123(6357))? = 125 gal/min > 120 gal/min required, ~. o.k. 


Another popular name for the self-venting gravity flow side outlet 
arrangement is overflow arrangement. Another popular name for the 
bottom outlet arrangement is drain arrangement. The entrance to an 


overflow arrangement is always above the horizontal surface. The 
entrance to a drain arrangement is always flush with a horizontal surface. 
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Horizontal Fully Vented Flow 


In many instances gravity flow is desirous in horizontal pipe lines running partially full. Obviously 
this cannot occur unless the line is slightly inclined to provide static head sufficient to overcome the 
frictional losses. This application is termed near-horizontal line flow. Hills? has recommended a 
minimum slope of 1:40 for these process situations with adequate free area remaining in the pipe to 
allow counter-current gas flow. 

While it has been stated that the theoretical maximum flow rate occurs in near-horizontal 
pipes at a flow depth of y = 0.9384, sufficient free area must be reserved in the upper region of the 
pipe to preclude gas carryover in the liquid and to allow for gas phase counterflow. Liquid depths in 
these two phase flows should not exceed more than '^ the inside pipe diameter for nominal pipe up 
to 8 inch NPS and not more than 74 of the diameter for larger pipes. 

A hydraulic jump, see Figure 8, will occur at the base of vertical runs when near-horizontal 
pipe lines flowing partially full turn to a vertical orientation, and then resume a near-horizontal 
orientation downstream. As long as there is sufficient distance to the first flow obstruction 
encountered in the near-horizontal pipe, the surge associated with this hydraulic jump is of no 
consequence. The jump normally occurs at approximately 10d downstream of the change in the 
flow orientation from vertical to near-horizontal. 


Pressurized Flow 


Introduction 


From confined space automotive applications to mammoth systems in industrial facilities, 
fluids are transferred by piping to and from and through the systems that are now taken for granted 
as part of our daily lives. Needless to say, a huge amount of energy is expended in accomplishing 
these transfers. Pressurized flow systems account for nearly 2096 of the world's electrical energy 
demand and range from 25-50% of the energy usage in many industrial plant operations." 


In the context of this course, pressurized flow means 
flow that is induced artificially by a pump, compressor, 


or a higher pressure source such as a boiler. 
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Design Approach versus Analysis 


Flow equations published in the open literature are invariably based on determining pressure loss 
due to fluid friction based on a known pipe size. As such, they represent an analysis approach. In 
many cases the requirements placed on a system are established 
in terms of a maximum allowable pressure drop or energy loss, a 
= D required volume or mass flow rate, the fluid properties, and the 
e material composition of the pipe to be employed. Accordingly, a 
suitable pipe size must be determined which will meet these 
requirements. As shown in previous sections of this course, situations where the determination of an 
optimum pipe size is required present a true design approach that can necessitate the rearrangement 
of the normal flow equations in terms of d. 


Estimation of Pipe Size for Liquids 


Adams" developed two mathematic equations which very closely describe the pressurized flow of 
liquids in typical pipes within the realm of normal supply pressures. They offer a simple way to do 
preliminary pipe sizing. If used for typical flow situations for normal fluids, they provide a quick 
size determination. 

The Adams equations in their original form are cubic in terms of pipe diameter, and as such 
are more suited for determining the approximate flow corresponding to a known nominal pipe size. 
However, with manipulation they can provide a quick means to determine an approximate nominal 
pipe size. A least squares geometric regression analysis of Adams’ equation for lower flow rates 
yields the following correlation equation: 


d = 0.250 [for O < 100 gal/min] [22] 


Adams’ equation for higher flows, rewritten in terms of d to allow for a direct solution yields, 


Q 1/3 
d= (2) -2 [for Q > 100 gal/min] [23] 
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Example 6 Quick Determination of Pipe's Flow Capacity 


Problem: Nominal pipe sizes have yet to be indicated on a preliminary P&ID (process and 
instrumentation diagram). Two separate process lines with non-viscous liquid flow rates of 85 gpm 
and 260 gpm are shown. What are the probable nominal pipe sizes for these two lines? 


Given: Q = 85 gal/min ; Q = 260 gal/min. 
Find: A target pipe diameter d for the stated flow rates. 


Solution: 

® Because handbook or table data are not available, use the reformations of the equations given by 
Adams" to make quick determinations. In the first instance the stated flow rate is less than 100 gpm 
-. use equation [22]: 


d = 02540 
d = 025485 
d - 23 inches 


^. A preliminary nominal pipe size for the first pipe line would be 275 inches. 


@ In the second instance the stated flow rate is greater than 100 gpm, .. use equation [23]: 


Q 1/3 
= |= 22 
4- (2 
| 260) 1/3 
EE =2 
g 1.2 
d = 4 inches 


-. A preliminary nominal pipe size for the second pipe line would be 4 inches. 


Pressure Conservation 


Invariably, in the case of many process fluids, the size determination becomes one that is driven by 
pressure conservation - that is to say, the size that will result in a specified minimum pressure value 
at some terminal point. System design problems usually contain too many unknowns to allow a 
direct solution. The resulting design approach can require a procedure known as iteration; a trial- 
and-error solution method in which a trial value is assumed for a given unknown, thus allowing the 
computation of a related unknown. The iteration procedure provides a means of checking the 
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accuracy of the assumed trial value and also indicates the new trail value to be used if an additional 
computational cycle is required. The following example illustrates the application of this design 
procedure. 


Instructors Note: While extremely useful, automated 
computational methods often preclude a clear understanding of 
the underlying engineering. The example that follows employs a 
classical method of design which was required prior to the 


availability of powerful desktop computing resources. It is 
presented here solely to illustrate the concepts of fluid flow 
relationships and their importance in selecting pipe size. 


Example 7 Optimum Pipe Size for Pressure Conservation (Classical Method) 
Problem: |n a petrochemical unit operation, para-xylene [C&H4(CH3);] at 86° F is required to be 
transferred at a rate of 100 gallons per minute through a schedule 40 steel pipe segment which has a 
length of 188 feet. The p-xylene has density and viscosity values of 53.6 Ib/ft? and 0.6 cP 
respectively. What is the smallest nominal pipe size that can be installed that will result in a 
pressure drop through the pipe segment not exceeding 2 psi? 


Given: Q — 100 gal/min (0.223 ft?/sec) ; L = 188 feet of steel pipe (e = 0.00015 feet) ; 
p = 53.6 lb/ft? ; u = 0.6 cP (o = 7.5 x 10° f/sec) ; AP = 2 Ib/in? (hy= 5.37 feet of p-xylene) 


Find: The smallest pipe diameter d for the stated process requirements. 


Methodology: The flow velocity V and the pipe diameter d are unknown so that neither the 
Reynolds number Ne. nor the relative roughness ¢/D can be calculated. The procedure will be to 
assume a value of fand calculate a corresponding value of d using the energy loss equation and the 
continuity equation. A Reynolds number and a relative roughness will then be calculated based on 
this value of d. Using the Moody diagram (similar to Figure 2), a new value of f will be obtained. 
This procedure will be repeated unit the value of fis repeated (or remains appreciably unchanged) 
and all of the flow equations are satisfied. During the design procedure, care must be taken to 
remember that classical flow equations employ the variable D (in feet). 
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Solution: 

© The flow velocity, Reynolds number, and the relative roughness are all dependent on the pipe 
diameter. Without these values the friction factor cannot be determined directly. Iteration is required 
to solve this system design problem since there are too many unknowns. 


@ Using a form of the continuity equation [V = Q/A] and the flow area formula [4 = xD"/4], 
rewrite the Darcy-Weisbach energy loss equation to appear as 


LE L 160° SLO f 
h= f Rz f 274 B 2 5 
g D z^D (2g) mg D 


Q Rewrite the resulting equation above in terms of D, insert the known values, and evaluate: 


 SfLQ'  (8)188(0223) 


D - 
m’ (3217)(537) 


3 f = 0049f 
gn h; 


® Use the continuity equation and the known quantities to write a form of the Reynolds number 
containing D: 


y -2 , 40 , (40295 _ 3.79 x10* 
F^ v  mvD 75xl0*zD D 


© Using the absolute roughness value for steel, the relative roughness in terms of D becomes 


E 0.00015 


D D 


© A trial value of f must now be assumed. Since both Nre and £/d are unknown, the Moody 
diagram yields no specific help with the assumption of the initial value. Referring to the left 
abscissa of Figure 2, assume a mid-range value of the Colebrook friction factor of 0.03. 


@ Compute the value of D from the Step 3 equation based on the initial assumption of f 
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D*= 0049f 
D = [(0049Y003)] ^ = 0.271 feet 


For this value of D, the following values can be computed 


4 4 0.00015 . 0.00015 
_ 379x10" _ 379x10 — | d ; Ee - = 0.00055 


Na, = D ~ 0271 D D 0271 


Q Figure 2 (or a more accurate version thereof) yields a Colebrook friction factor of f= 0.02 
corresponding to the above calculated Reynolds number and relative roughness. Repeating the 
process for f= 0.02 

D = [(0.049)(0.02)]"" = 025 feet 


| 379x10 _ T 
R^ ^25 —^* 
€ 000015 _ 35956 
D 02 |^ 


The Moody diagram yields a Colebrook friction factor of f= 0.019 corresponding to the above 
calculated Reynolds number and relative roughness. Since this value of f'is not appreciably different 
from the previous iteration, D — 0.25 feet (d — 3 inches) is probably the optimum size. From Table 
13, select as the optimum pipe size a 3 inch NPS Schedule 40 pipe with d= 3.068 inches . Asa 
check, the actual pressure drop for this selected pipe size will be computed to insure that the 
maximum allowable pressure drop has not been exceeded: 


0.223)(4 
= a = 435 ft/sec and 
z (3068/12) 
2 2 
2 fee soe) - 4.21 feet p-xylene = 1.57 Ib/in’ .. o.k. 


"=F Dag” 975 G068122)3217) 
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Pipe Size via Fluid Velocity 


The most common pipe size selection criteria is most usually based on fluid velocity. Years of 
experience have generated rules-of-thumb regarding normal velocities for given fluids. Knowing 
the flow rate, optimum pipe size is back-calculated from the traditional continuity equation which 
can be rewritten in the form 


40 


47 Wav 


[24] 


The generally accepted body of data for normal velocities is presented in the Table 8. 


TABLE 8 Widely Published Recommended Flow Velocities 


Mean Flow Velocity 
( feet per ( feet per 


Service Minute) second) 
Boiler Feed Water 500 — 600 8-10 
General Water Service 200 — 500 3-8 
Process Water 300 — 600 5-10 
Pump Discharge 300 — 600 5-10 
Pump Suction 200 — 500 3- 8 
Saturated Steam (high pressure) 6,000 — 10,000 100 - 165 
Saturated Steam (low pressure) 4,000 - 6,000 65 - 100 


Superheated Steam (high pressure) 10,000 — 15,000 165 - 250 


For low-viscosity liquids economic optimum velocity is 
typically in the range of 6 to 8 feet per second. For gases with 


density ranging from 0.015 to 1.25 lb/ft? the economic 
optimum velocitv is about 130 down to 30 feet per second. 


The term optimum velocity has often been incorrectly limited to mean the velocity for a 
certain flow rate and pipe size (flow area) which results in an economic balance between hardware 
and energy costs. Optimum velocity may also very well mean the velocity which will maintain a 
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certain degree of agitation; or will produce an efficient or improved heat transfer; or may limit 
abrasion or conserve the suspension of solids. 


Pipe Size based on Fluid Velocity Effects 


Not all fluid particles travel at the same velocity within a pipe. The shape of the velocity curve (the 
velocity profile across any given section of the pipe) depends upon whether the flow is laminar or 
turbulent. If the flow in a pipe is laminar, the velocity distribution at a cross section will be 
parabolic in shape with the maximum velocity at the center being about twice the average velocity 
in the pipe. In turbulent flow, a fairly flat velocity distribution exists across the section of pipe, with 
the result that the entire fluid flows at a given single value. Figure 12 helps illustrate the above 
ideas. The velocity of the fluid in contact with the pipe wall is essentially zero and increases with 
the distance away from the wall. 


Laminar Flow 
N,.<2,000 
Re 


Laminar Flow Turbulent Flow 


Figure 12 — Velocity Profiles (adapted from Reference 13) 


In many process applications pipe size is dictated by fluid velocity effects. Fluid velocity 
can affect noise generation, vibration, erosion, and the fluid's transport efficacy. Flow velocity is 
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normally maintained high enough such that liquid flows are swept along avoiding the slug flow 
pattern. 

This section of the course addresses two important considerations of fluid velocity and 
therefore pipe size selection. These considerations are flow induced erosion and settling velocity. 

For a given required mass flow rate, pipe size can have a significant bearing on internal pipe 
erosion. An important consideration in the selection of pipe to provide pressure resisting strength is 
the tolerable amount of erosion, in terms of loss of pipe wall thickness. Erosion occurs in both 
single and multiphase flow patterns. With regard to single phase flow, studies have shown that gas 
flows are generally more damaging than liquid flows.'* With that said, flow of concentrated liquid 
acids have been found to be extremely erosive. 

It logically follows that erosion in straight pipes is less severe than that in elbows and other 
fittings for a comparable set of operating conditions. In elbows fluid can possess enough momentum 
to traverse the flow streamlines and make direct impinge against the pipe wall. Erosion can also 
occur in straight sections of pipe even though there is no mean velocity component directed toward 
the wall. This occurs when turbulent fluctuations in the flow provide the matter with momentum in 
the radial direction, forcing it into the pipe wall. 

In the oil and gas industries, the presence of contaminants in produced oil and natural gas 
represents a major concern because of the associated erosive wear occurring in various flow 
passages. The commonly used practice for controlling erosion in gas and oil producing wells is to 
select a pipe size to limit production velocities following the provisions of American Petroleum 
Institute. Recommended Practice 14E. This document contains a simplified, very conservative 
formula to calculate a threshold velocity, and thereby a corresponding pipe size, below which an 
allowable amount of erosion occurs. Rewritten in terms of d with gas flow Q in ft/sec this 
expression becomes, 


[25] 


While providing the advantage of computational ease, this approach has some 
disadvantages. The variable c is an empirical constant, intended to represent the normally occurring 
corrosive and erosive effects of entrained contaminants in the fluid. 
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A misconception of the API erosive velocity formula is 
that values of c account for varying degrees of sandy 


petroleum flows. In fact, assumption of values for c is 
intended to allow for a range of clean service conditions. 


The original formula from which Equation [25] is derived has been found to be very 
conservative for relatively clean fluids and inapplicable for corrosive/sand containing processes. 
The definition of its use and the validity of the assigned values of 100, 125, 150, and 200 for c have 
generated a host of supplemental third party guidance literature." The value of c is irrespective of 
volume or weight fraction of gas-to-liquid in the petroleum. Multiphase flow is accounted for in the 
formula by the incorporation of the variable representing the density of the fluid mixture. In real- 
world applications many factors, in addition to the density of the mixture, influence the erosion rate. 
These include flow geometry, particle size, and Reynolds number. Shirazi et al. have presented a 
method to calculate threshold velocities to overcome these limitations by accounting for many of 
the physical variables in the flow and erosion processes and includes a way to predict the maximum 
penetration (pipe wall thickness reduction) rate for sand erosion. In reality erosive material can 
travel in a non-uniform, transient slug flow pattern, which settles and forms deposits. As the 
deposits increase, the flow cross-sectional area decreases, and the resulting erosive velocity 
increases, entraining and slugging through the deposits that have previously settled. 

Gas erosional velocity and pipe size selection will be discussed in more detail later in the 
course. 


Solids Suspension, Settling, and Sedimentation 


If a fluid, such as water, is flowing, it can carry suspended particles. The settling velocity is 
the minimum velocity a flow must have in order to transport, rather than deposit, sediments, and is 
given by Stokes Law. 

If the flow velocity, often in these instances referred to as shear velocity, is greater than the 
settling velocity, sediment will be transported downstream as suspended load. As there can be a 
range of different particle sizes in the flow, some will have sufficiently large diameters that they 
settle on the lower pipe wall, but still be displaced downstream. In many cases of process design, 
the terminal velocity of spherical particles (the settling velocity) must be ascertained. In most cases 
the application of Stokes’ law is appropriate; however, in some cases other flow regimes apply 
which require a different form of the equation to determine settling velocity. 
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McCabe and Smith“ give a method for determining the correct flow regime using physical 
data. A dimensionless constant K is employed which has the value of: 


[26] 


e, (Pp - eo) 
u 


K= ZETA ; 


Flow regimes are categorized in ranges of either Stokes range, intermediate range, or 
Newtonian range. 


TABLE 9 Constants b and n in Equation [27] 
for three operative flow regimes 


Law Range b n 
Stokes K<3.3 24.0 1.0 
Intermediate 3.3 > K < 43.6 18.5 0.6 
Newton’s K > 43.6 0.44 0 


Caution: computed values of K exceeding 2,360 
result in a regime which exhibits abrupt changes in À 


the drag coefficient with minimum changes in 
velocity, nullifying the equations presented here. 


After the appropriate flow regime has been determined, the sediment or solids single particle 
settling velocity can be found from 


lin len 1/(2-n) 
Pa D, (Pp Pr) 4387(1/12) [27] 
SOO |. bu” pt (6.72x104)” 
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using the constants b and n given in Table 9. Single particle settling velocity by itself is of limited 
usefulness since practical applications generally involve agglomerations. Fortunately, Perry! has 
provided a relationship which accounts for hindered settling due to the interaction of proximate 
particles 


Vy =  Vgü- o)" [28] 


Where m has the value shown in Table 10: 


TABLE 10 Single Particle Settling Velocity Modification 
Exponent for Hindered Settling Equation [28] 


Reynolds No. Range m 
Nre <0.5 4.65 
0.5 2 Nre < 1300 4.375 Nre "5" 
Nre > 1300 2.33 


For a rigorous treatment of this topic it is recommended that Reference 18 be read. Barring 
its availability, the practical concept of solids suspension conservation is best illustrated by the 
example that follows. 


In fluid mechanics two-phase flow occurs in a system 
containing gas and liquid. The term multiphase flow tends to 


be used to characterize two-phase flow in cases where the 
two phases are not chemicallv related (e.g. dustv gases). 
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Example 8 Conservation of Suspended Solids in a Transfer Pipe 


Problem: lt is proposed to transport iron particulate through a pipe via a gas stream at a rate of 

50 ft?/min. The iron particles, which will constitute a mean volume fraction of 10% in the resulting 
mixture, have been determined to be of uniform size and fairly spherical, measuring 0.01 inches in 
diameter. The particle density is 500 Ib/ft?. The transport fluid (gas) density and viscosity are 0.08 
Ib/f? and 0.02 cP respectively. Without regard to pressure drop or erosion, what is the probable 
nominal diameter pipe that can be used to accomplish the multiphase transport of this metal? 


Given: Q = 50 ft?/min ; c = 0.10 ; Dp = 0.01 inches ; pp = 500 Ib/f? ; pr = 0.08 Ib/fP ; u = 0.02 cP 


Find: A target pipe diameter d for the stated process application. 


Solution: 

© In order to have unsedimented transport, the shear (transport) flow velocity must be greater than 
the suspended iron particle hindered settling velocity. First determine the appropriate flow regime 
by computing the McCabe-Smith constant using Equation [26], 


E 1/3 
K= usi, reer - 48100] 


(0.08)(500- 0.08) | '? 
0.02? dn 


Since K = 16, from Table 9 it can be seen that the intermediate flow regime is applicable and 
consequently the constants b and n are assigned the values of 18.5 and 0.6 respectively. 


@ Next, calculate the single particle settling velocity from Equation [27], 


n 1/(2=n) " 
y.- n? (P, = Pr) 43.87 (1/712)! E (aon (500- 0.08) 43.87 (ine 1/(2—-0.6) 
: b u^ p? (672 x10^)" ~ | (185)(0.02)^5(0.08)- 95 ^ (6.72 x10") 
Vs = 12 ft/sec 


© Compute the particle Reynolds number using 


T 124 D, Vs Pr 7 (124)(0.01)(12)(0.08) 
Re 7 E 
u 0.02 


= 60 


(the coefficient of 124 converts inconsistent units) 
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€) Since Nre = 60, from Table 10 find m from the expression 


m = 4.375 Nz’ = (437560) 9995 = 3.06 


© Determine the hindered settling velocity of the iron solids using Equation [28] 
Vor= V;ü-c)- 120-010)95- 87 ft/sec 
© Select a suitable pipe diameter which will produce a flow velocity that exceeds V5; . Use the re- 


arranged continuity equation previously presented to arrive at a target size. Be careful to convert 
quantities to insure consistent units. 


d- 40 - [ee ecm) - ene 
Y UV sry (min)(z)(8.7ft)(60sec)(ft^ ) 


@ Since this inside diameter will result in a flow velocity equal to the hindered settling velocity, a 


value actually smaller than 4.2 inches must be used in order to produce a flow velocity that exceeds 
the hindered settling velocity. 


Compressible Flow 

Flow where the density varies by more than 5% is considered compressible. An incompressible 
flow is flow in which the density of the fluid is constant or nearly constant. Liquid flows are 
normally treated as incompressible. Unlike liquid steady-state flow where the mean flow velocity 
remains relatively constant for a given pipe size, gas flow velocity can vary along the pipe’s length. 
This is due to compressibility effects, and variation in pressure and temperature. 


The Compressibility factor Z is a dimensionless parameter less than 1.00 that 
represents the deviation of a real gas from an ideal gas. At low pressures and 
temperatures Z is nearly equal to 1.00 whereas at higher pressures and 


temperatures it may range between 0.75 and 0.90. The actual value of Z at 
any temperature and pressure must be calculated taking into account the 
composition of the gas and its critical temperature and pressure. 
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In a given pipe size of fixed cross sectional area were the mass flow rate of the gas is 
constant, pressure, temperature, and density changes in the gas will modulate the velocity. If the gas 
flow at some upstream point in the pipe is at a velocity less than Mach 1 (subsonic), the pressure, 
temperature, and density of the gas can decrease to a point were the velocity becomes Mach 1. It is 
at this point within the pipe where the flow becomes sonic and any intended means of upstream 
flow control, e.g. pressure relief valve, is rendered ineffectual. This flow condition is said to be 
choked, a term that it is semantically misleading because the flow with this condition is far from 
being totally obstructed. 


CHOKED FLOW POINT 


FLOW RATE 


1 


Figure 13 — Sonic Flow 


Compressible flow problems can be complex and very difficult to model. For simplification, 
often times if the anticipated compressible flow pressure and/or temperature changes are small, the 
flow is analyzed as incompressible. Additionally, transmission of gas in overland pipe lines is 
generally considered to be carried-out isothermally. In contrast, shorter length compressible gas 
flows are assumed to occur adiabatically (no heat transfer). Of course in reality flow conditions 
probably fall somewhere between these two thermodynamic processes. 

Owing in great part to the natural gas industry, there is an extensive number of empirical 
flow equations which are in place which attempt to model compressible flow in pipe. The general 
gas flow equation rewritten in terms of d is, 
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i 70.4 
29 
" Q [29] 
P- P 
8306. cary f 


Where the subscripts 1 and 2 denote upstream and 
downstream conditions, respectively. 


Static head of gas column (+ or — ) 
correction factors must be introduced to 
the general gas flow formula when 


significant topographical elevation 
differences exist between the source (inlet) 
and terminus of an overland gas pipe line. 


Important considerations when dealing with gas flow formulae to determine optimum pipe 
size are: 


Pressures and temperatures are always taken at their absolute value; 
Volumetric flow rate is in standard ft*/day, often times abbreviated SCFD; 
Gas specific gravity is relative to air at standard conditions; 


ee a eae 


Pipe lengths are in miles in lieu of the customary feet. 


According to Menon? the velocity in an isothermal compressible gas flow at any point is 
described by, 


ZTQ 

Feoxlü" = | 30 

T pg EN 
Where Q is in standard ft*/day 
and the limiting erosional velocity for such a flow is given by the empirical formula, 
ZTR 
-] — 31 

Me My 29PG BY 
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A The variable R in Equation [31] is the universal gas constant, not the hydraulic radius. 


Equating expressions [30] and [31] and evaluating the constants results in the following useful 
mathematical relationship to select a pipe size which mitigates gas flow pipe wall erosion: 


ZTG 


d = ooo Jo | 22) B [32] 


Example 9 Optimum Pipe Size for Gas Flow with Limiting Velocity 


Problem: Natural gas is to be transferred at the rate of 3.8 x 10° standard cubic feet per day through 
11 miles of transmission pipe line. Assuming the average pressure is 900 Ib/in? and the flow is 
carried-out isothermally at 70? F, what pipe size should be employed to minimize erosion to the 
pipe wall? The compressibility factor and specific gravity for this natural gas have been determined 
to be 0.85 and 0.60 respectively. Calculate the maximum allowable erosional velocity for these flow 
conditions and estimate the actual fluid velocity for the selected pipe size. 


Given: Q =3.8 x 10? f/day ; L = 11 miles ; P = 900 Ib/in? ; T= 70° F ; Z=0.85 ; G = 0.60. 
Find: The smallest nominal pipe diameter d for the stated duty which will minimize erosion. 


Solution: 


(D Convert pressure and temperature to absolute values: 
P =900 + 14.7 = 914.7 Ib/in absolute 
T= 70 + 460 = 530° R 


Q Use the given values and Equation [32] to find that, 


Zug 0.85)(530)(0.60) | ^ 
d = ooo Jo | 22) - vom ssa! | S30 i = 14.373 inches 


Q From Table 13 on page 77 the closest standard weight nominal pipe size is a: 
16 inch NPS schedule 30 with d = 15.250 inches 


® The limiting erosional velocity is found from Equation [31]: 
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ECT {Cocos — 
a 29PG ` (299147 (0.60) ^ ^" em 
®© The resulting velocity for the selected optimum pipe size is, 
ZT 0.85)(530)( 3.8.x 10° 
V = 6x10” 9 = 6x10. N30) > ). 48.3 ft/sec 
Pd (914.7)(15.250) 


© Since V « V, (48.3 < 55.1), then the selected pipe size is satisfactory. 


The Weymouth gas flow equation is the oldest (1912) and most common of the empirical 
formulas. A variant of its original form, rearranged to allow a direct determination of pipe size is, 


0.1876 


2 
] | WGL [3] 


7868(P? - P?) 


Example 10 Gas Pipe Line Size based on the empirical Weymouth Formula 


Problem: Gas is to flow at the rate of 85,000 Ib/hr through an overland transmission line. The 
guaranteed delivery pressure 9/4 miles downstream is 20 Ib/in?. The compressor will discharge the 
gas, which has a specific gravity of 0.80, at 200 Ib/in?. It has been estimated that the average 
operating temperature will be 60? F. What size pipe should be used for this service? 


Given: W = 85,000 lb/hr ; L =9.5 miles ; P; = 200 lb/in’ ; P = 20 Ib/in? ; T= 60? F. 


Find: A suitable pipe diameter d for the stated utility application. 


Solution: 

(D Convert pressures and temperature to absolute values: 
P, =200 + 14.7 = 214.7 Ib/in? absolute 
P; =20 + 14.7 = 34.7 lb/inņ? absolute 
T= 60 + 460 = 520° R 


Q Use the given values and Equation [33] to find that, 
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g | WGL 


° | (85,000)? (080)(9.5) 
7868( P? - P2) 


0.1876 
= = 3.970 inches 
7868(214.7? - = 


Q From Table 13 on page 77 selecta 4 inch NPS schedule 40 pipe with d = 4.026 inches 


Two-Phase Flow 


The majority of two-phase (liquid-gas) flows manifest in patterns known as bubble, plug, slug, 
wavy, and dispersed. These are shown in Figures 4 and 5. Multiple phase flows, even when 
restricted to simple pipeline geometry, are in general quite complex, and several features may be 
identified which make them more complicated than single-phase flow. As shown in the two figures, 
because of phase density differences, vertical flow patterns are different from horizontal flow 
patterns, and horizontal flows are generally asymmetric. Even when phase equilibrium is achieved 
by good mixing in two-phase flow, the changing equilibrium state as the pressure drops with 
distance, or as heat is added or lost, may require that inter-phase mass transfer, and changes in the 
relative amounts of the phases, be considered. 

Because two-phase flow applications and their flow equations can be so challenging with 
regard to selecting pipe size directly, an indirect approach is often employed in which a trial pipe 
size is selected and the corresponding pressure drop is evaluated. This cycle is repeated until a size 
is found that results in an acceptable energy loss for the required flow rate. Lockhart & Martinelli” 
developed a semi-empirical method that conservatively approximates a two phase pressure drop that 
is applicable to the four major two phase flow patterns of stratified, wavy, slug, and dispersed. 
Although still an excellent approximation, the method has been found to under estimate the actual 
pressure drop for the annular flow pattern. The Lockhart & Martinelli correlation is based on the 
application of a factor derived from the independent single phase pressure drops of the two phases. 
This factor is applied to the single phase pressure drop of either phase to arrive at the mixture 
pressure drop. Like most other analytical exercises, determination of an appropriate friction factor 
presents the single biggest burden. 

The Lockhart & Martinelli method uses normal flow equations to compute common values 
associated with single phase pressure drop as if each phase was flowing singly in the pipe. The 
Blackwell?'/Lockhart & Martinelli method then applies Equations [34], [35], [36], and [37] to arrive 
at the total mixture pressure drop. The subscripts / and g refer to liquid and gas respectively: 
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+ 
sl NE JP [34] 


Two-ph ixt locity”: 
wo-phase mixture velocity V = 0.0509 d? 


Lockhart & Martinelli two-phase flow modulus”: X = AP [35] 
g 
2 2 2 
Blackwell correlation of Y, from Ref. 22: y = [oem *A,lnX ) [36] 


Where, Ao = 1.4659 ; A; = 0.4914 ; 4; = 0.0489 ; A; = -3.487 x 10* 


Lockhart & Martinelli two-phase mixture pressure drop: AP = AP, Y, [37] 


A The above equations are applicable to the fully turbulent flow regime. 


Rearrangement of Equation [34] allows for the computation of a suitable pipe diameter when a 
mean fluid mixture velocity is known: 


[38] 
d = 0.2256 


The best way to illustrate the Blackwell/Lockhart & Martinelli methodology is through a worked- 
out example: 


Example 11 Optimum Pipe Size for Two-Phase Flow with Specified Velocity 


Problem: A fluid mixture of gas and liquid must be conveyed at a mass flow rate of 650,000 Ib/hr. 
The weight fraction of the liquid which has a density and viscosity of 33.5 Ib/f? and 0.10 cP 
respectively, is 46%. The gas has a density of 2 lb/ft? and a viscosity of 0.01 cP. If a minimum fluid 
mixture velocity of 66 feet per second is required to preclude damaging slug flow, approximate the 
probable diameter of pipe that must be installed to meet these conditions. Calculate and evaluate the 
resulting pressure drop assuming this steel pipe 1s 100 feet long. 
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Given: W = 650,000 lb/hr ; W, = (0.46)(650,000) = 299,000 lb/hr ; .. W, = 351,000 lb/hr ; 
pi = 33.5 Ib/f ; 14 = 0.10 cP ; p, = 2.0 lb/ft? ; uz — 0.01 cP; V = 66 ft/sec ; L = 100 feet of steel 
pipe (e = 0.00015 feet). 


Find: The target pipe diameter d and a corresponding pressure drop for the given process flow. 


Solution: 


O Using Equation [38] and the given quantities, calculate an inside pipe diameter to be, 


AM i4 351,0007 . 299,000 
Pg /Pi 2t 33.5 
e EET jaa 
V 66 


d = 0.2256 = 11.93 inches 


@ Knowing that D = d/12 = 11.93/12 =0.994 feet, calculate the relative roughness e /D to be 
0.000151. 


®© Assume the gas is flowing alone and calculate a Reynolds number? of, 


351,000 
(11.93)(0.01) 


e 


W 7 
Mac 63177 = 631 - 186x10 


@ Knowing these values of ¢ /D and Nre , use a Moody chart, internet applet, or any accepted 
empirical formula to compute a Colebrook friction factor of f= 0.013. 


® Calculate the gas pressure drop using a form? of the Darcy-Weisbach formula, 


LW 0.013)(100)(351,000)? 
£ - 336x190 5f Ee, 


AP, = 336x106 ZE : 
$ d (11.93) (2) 


Pg 


- 11134 Ib/in? 


© Now assume the liquid is flowing alone and repeat Steps ©, ©, and © and calculate: 
Nre = 1.59 x 106 
f = 0.014 
AP, = 0.0515 Ib/i? 


( Next, find the Lockhart & Martinelli two-phase flow modulus using Equation [35]: 
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3n AP, S — 
7 ~ V11134 ` 


Using this value of X, calculate Perry's flow parameter Y, from Equation [36], 


Y= 


2 
: oes | = 5237 


GO Now compute the mixture pressure drop using Equation [37] and evaluate its suitability: 


AP= AP, Y; = (1.1134)(5.237) = 5.83 Ib/in? 


Pressure Relief Valve Discharge Pipes 
One of the most important pipe size selection exercises is that associated with pressure relief valves. 


The inlet, and especially the outlet piping of pressure relief valves, is critical to their intended 
proper performance. The discharging fluid friction in the discharge piping produces back pressure 
on the relief valve; the larger the pipe diameter, the lower the resulting back pressure. 


Relief valves are designed to work with back pressure from 10% 
(conventional valves) to 40% or more (balanced bellows valves) of 


their set pressure. For example, a conventional valve relieving at 
150 psig will work reliably if the back pressure is limited to 15 psig. 


The simultaneous discharge from multiple relief valves makes line sizing more complex. 
Relief discharge pipes sized for protecting tanks against external fires requires the assumption that 
all valves associated with a manifold will relieve at once.™ Safety relief systems where the relieving 
conditions are based on other criteria, e.g. runaway chemical reactions, must be analyzed to 
determine how many devices could conceivably relieve simultaneously. The actual task of pipe 
sizing can begin once it has been determined which valves will relieve simultaneously. 
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The relief valve sizing procedure often times necessitates the 
selection of a relief valve flow orifice which results in a larger 
relieving flow rate than that required by the process. Pipe sizing 
must be carried out for the flow rate corresponding to the actual 


relief valve orifice size selected at set pressure, and not at the 
originally envisioned required process relieving rate. See 
PDHcenter.com course no. M112 for a detailed treatment of the 
selection and sizing of pressure relief valves. 


Two-phase flow through pressure relief valves is common-place. Assumption of a complete 
vapor condition will result in conservative results where multiphase flow is apparent and the total 
flow rate is known. Going further, assumption of isothermal flow 
conditions has shown to result in higher pressure drops (more 
conservative estimates) than an adiabatic process assumption. Actual 
flow more than likely lies between these two thermodynamic 
processes. Because two-phase pressure drop prediction models are 
inherently uncertain, the conservatism is welcomed. 

The rapid expansion of a single phase liquid through an 
actuated pressure relief valve can, and often does, result in flashing 
flow in the discharge pipe. Contrasted to two-phase flow which 
denotes steady-state conditions, flashing is the term used to describe 
the special category of liquid vaporization which involves a fluid’s 
rapid phase change from liquid to vapor without the return of the 


fluid to the liquid phase. 

Relief valve discharge manifolds are used over separate 
dedicated pipe runs to minimize costs and protect the environment. Pipe size section in relief 
manifolds is critical to the reliable operation of a safety system. While time consuming, manual pipe 
sizing is straightforward. Each pipe segment in the manifold network is analyzed based on the flow 
rates summed from contributing upstream relief devices. It is often necessary to perform the 
computations using different scenarios. The proper way to size the piping is to proceed in reverse 
(upstream) from the point where the manifold discharges to a collection element or atmosphere. The 
best approach is to make estimates for trial pipe sizes based on the equations presented in this 
course. Detailed calculations are then performed to determine the pressure drop through each pipe 
segment. The allowable back pressure for each relief valve is compared with the back pressure 
calculated at the valve. Technical reasoning is used to adjust pipe sizes, larger or smaller, in the 
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network; the calculations are repeated and the procedure continued until all back pressures are 
satisfactory. Keep in mind that critical (sonic) velocity must not be exceeded in a given pipe 
segment in order to preclude choking flow. 

As in single pipe lines, actual flow in relief valve discharge manifolds also normally takes 
place somewhere between adiabatic and isothermal conditions. The assumption of an isothermal 
model is a conservative approach in that isothermal equations yield higher manifold backpressures 
for both subsonic and sonic flow.” 

Equation [43] presented on page 60 offers an iterative, manual approach to selection of pipe 
sizes for flare and relief discharge manifolds. Students who are less interested in the development of 
the methodology may wish to proceed there now. 

Kandell5 through the combination of sonic and actual gas flow equations and the 
incorporation of the ideal gas law showed that the Mach number for isothermal flow is given by: 


W [ZT 
M = 1336x102 YT EUR [39] 
Ww 


Knowing that A = zd’/4 allows the transformation of Equation [39] to appear as: 


M = 000244951 —— E 40 
m Pd? \ M, mié. 


which permits the direct insertion of an assumed or selected pipe size. Rearrangement of Equation 
[40] in terms of d allows for the direct selection of a pipe size with an assumed upper limit for the 
Mach number: 


0.5 
WJZT|M,, 
TA [41] 
408.245 MP 
[with the value of M usually assigned in the range of 0.6 to 0.8] 
For isothermal flow is can be shown that. M A= M; P, [42] 


An adaptation from Kandell and Lapple?' shows that for isothermal flow based on upstream 
conditions: 
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P 2fL P 
1 -= m| rt + Ine +] [43] 


Where the subscripts 1 and 2 denote upstream and downstream conditions, respectively, for a given 


pipe segment. 


Relief discharge manifold networks contain too many unknowns to allow for a direct 


solution for pipe size selection. A rigorous design approach requires the procedure known as 


iteration which has been previously explained. The use of Equation [43] provides an iterative design 


approach and pipe size selection method for two-phase relief discharge manifolds. This is best 


illustrated by a worked out example which will follow momentarily. Here is the step-by-step 


analytical approach to the selection of pipe sizes for flare and relief discharge manifolds where 


compressible flow is evident and simultaneous valve actuation is probable: 


I. 


Construct an arrangement drawing of the discharge piping arrangement and assign 
location labels, A;B;C; etc. to each relief valve discharge point. Label the pipe 
intersections as well; 

Designate on the drawing the known relief flow rates, temperatures, molecular 
weights, and viscosities at the relief valve discharge points. Creation of a separate 
data table is also helpful; 

Based on the style of relief valve and its set pressure, compute the maximum 
allowable back pressure (MABP) at the discharge point of each relief valve. Add this 
data to the arrangement drawing and the data table; 

Calculate the combined flow rates downstream of each pipe intersection and add this 
data to the drawing and the data table; 

Estimate the fluid mixture properties downstream of each pipe intersection by the 
flow rate weighted average method using the following formulas: 


"2537 Yu) [44] , re Yu Y [45], 
i-l i=1 i-l i=1 

Hs Y su Ji.) Y «i. [46] 
i=l i=1 


Where n represents the number of constituent fluid streams. Add these fluid mixture 
properties to the drawing and the data table. Keep in mind that compressible flow 
equation variables must be evaluated with absolute values; 

Begin the analysis at the terminal point (or flare tip) of the manifold, and progress 
backward. In order to conduct the analysis, a system terminal point pressure must be 
assumed. This assumption is simplified if the terminal point of the manifold is 
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10. 


11. 


12. 


13, 


14. 


Is. 


16. 


vented to atmosphere at sea level: the value of the terminal pressure is 14.7 Ib/in? 
absolute. 

Based on the known (or otherwise determined) terminal point pressure and the 
estimated fluid properties in the final pipe segment, select a trial pipe size for this 
segment based on Equation [41]. Use good engineering judgment in the assignment 
of a value for the Mach number which is appropriate for the process. Keep in mind 
that it is desirous to limit fluid velocities such that M « 1 (subsonic) in order to 
preclude excessive vibration and unacceptable noise generation. In order to 
introduce a factor of safety, limiting the flow velocity at any point in the manifold 
corresponding to a maximum allowable Mach 0.6 is advisable. Design values range 
from 0.6 to 0.8; 

Using Equation [40] calculate the terminal point Mach number corresponding to the 
trial pipe size selected. If it is greater than 0.8, then the velocity is too close to sonic 
and the assumed pipe size is too small; 

Based on the selected trial pipe size, compute the pipe's relative roughness (c/D) and 
the pipe segment Reynolds number using an appropriate equation. In combination 
these values allow for the determination of the pipe segment's fluid friction factor; 
Substitute known values into Equation [43] to establish an iterative test equation 
which will allow for the solution of the upstream or inlet pressure for this pipe 
segment; 

Once the inlet pressure is computed, calculate the pipe segment's inlet Mach number 
using Equation [42]; if it is greater than 0.8, then the velocity is too close to sonic 
and the assumed pipe size is too small. Select the next larger nominal pipe size and 
repeat Steps 9, 10, and 11 until an acceptable Mach number is obtained; 

Working backward, the inlet pressure established in Step 10 should now be thought 
of as the outlet or downstream pressure for the preceding upstream pipe segment 
which will be analyzed. Similarly, the inlet Mach number determined in Step 11 now 
becomes the Mach number at the downstream point of the preceding pipe segment. 
Repeat this exercise beginning with Step 9, for each pipe segment in the manifold, 
working backward toward each relief valve, until the pressure at each relief valve has 
been established; 

The resulting calculated pressure at each relief valve is the back pressure seen by the 
relief valve. It should be compared to the MABP for that valve. The back pressure 
must be less than the maximum allowable; 

If the back pressure exceeds the MABP, a larger pipe size must be selected for the 
pipe segment under consideration. 

If there is a significant positive difference between the calculated back pressure and 
the MABP, the longest pipe segment should be decreased in size until the calculated 
back pressure approaches the MABP. 


The following example illustrates the application of the manual design procedure just 


outlined. As an alternative to this manual method, commercially available computer programs 
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which utilize looping iteration features to converge on solutions to implicit mathematical 
relationships handle complex flow problems like this quickly and accurately. 


Instructors Note: While extremely useful, automated 
computational methods often preclude a clear understanding of 
the underlying engineering. The example that follows employs a 
classical method of design which was required prior to the 


availability of powerful desktop computing resources. It is 
presented here solely to illustrate the concepts of fluid flow 
relationships and their importance in selecting pipe size. 


Example 12 Selection of Pipe Sizes for a Pressure Relief Discharge Manifold 


Problem: To minimize piping installation costs, the pressure relief discharge manifold shown in 
Figure 14a which ultimately vents out-of-doors, will be constructed in a chemical process plant. The 
worst case scenario for the process has been determined to be loss of cooling water; this will result 
in the simultaneous actuation of all four relief valves. Flashing compressible flow discharge 
conditions are documented. The known fluid properties and relief rates at each valve discharge 
point are indicated in the figure. Assuming fully turbulent isothermal flow conditions of ideal gases, 
what size pipe should be used for each segment of the manifold in order to insure that maximum 
allowable back pressures are not exceeded and subsonic flow conditions are maintained? 


Given: The process information shown in Figure 14a and assumption of fully turbulent isothermal 
relief discharge flow conditions with Z =1. Material: alloy steel (e = 0.00015 feet). Segmental 
lengths given are equivalent pipe lengths that account for fittings and entry and exit losses. 


Find: The smallest pipe size for each of the indicated pipe segments which will result in acceptable 
relief valve back pressures and universal subsonic flow. 


Methodology: The same analytical sequence presented in the course will be used in this example; 
however, the numbered steps in the solution that follows do not coincide with the step-wise 
enumeration in the course. In the interest of brevity, sample calculations are presented in lieu of 
repetitive computation details. Nevertheless, the final completed computational data and design 
results for the entire system are indicated in Figure 14b and summarized in Table 12. 
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o— [0 


L 7 1,118 ft 
CO,* HO L-100ft | CH 
W= 110,000 Ib/hr e W = 125,000 Ib/hr 
T=610°R e va T = 640°R 
My = 40 A M, = 84 
u = 0.0100 cP a=? a u = 0.0110 cP 
i 60 © n Of 
An (c) H d=? -M 
HT E— ` ET i 
È d=? L=125ft HE 
PE: [| L-300 ft T 
d=? o 
= 
= te) 
e = 
y he 
Ha e 
d Hu ; 
0 2 
$T |< 
i3 SO;,* H;O 
ao W = 75,000 Ib/hr 
T=580°R A CH; 
My = 60 W = 60,000 Ib/hr 
u = 0.0098 cP [ T = 800°R 


Tr M, = 54 
(F) u = 0.0130 cP 
Bal. Bellows PRV 


Set Press = 78 


Figure 14a — Relief Discharge Manifold Diagram for Example 12 
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Solution: 


© In relief valve parlance the maximum allowable back pressure (MABP) is equal to a fractional 
reduction of the valve's set pressure (generally 10% for conventional valves, 40% for balanced 
bellows valves). Sample Calculation: For the conventional valve at point E: 

MABP = (310)(0.10)+ 14.7 = 45.71b/in^ absolute 
Refer to Figure 14b for the calculated values of MABP for valves at points F, G, and H. 
@ Sample Calculation: By inspection, the mass flow rate downstream of point C will be: 

W = 75,0004 110,000 = 185,000 lb/hr 

Refer to Figure 14b for the calculated values of mass flow rate W, downstream of points B and D. 


®© Sample Calculation: The estimated fluid properties downstream of point C are: 


z r ` 75,000 + 110,000 « 46 
Equation [44]: M,, = ^n » (m M, ), ~ 750006 10000 t 
i=l i=l 60 * 40 


Y Y 75,000)(580) + (110,000)(610) 
: . T Š W, D W, " ( > , " o R 
Equation [45] A a 75,000+ 110,000 598 


By inspection, the mixture fluid viscosity will be u = 0.0099 cP. 


Refer to Figure 14b for the calculated values of fluid mixture properties downstream of points B 
and D. 


® By definition, the pressure at point A is P, = 14.7 Ib/in? absolute. Determine a trial pipe size for 
pipe segment AB based on P, and Mach 0.6 using Equation [41]: 


w |zT/M, |  [370000.(645)55.8] 
- EA - Bur = 18.69 inches 


408.245M P, (408.245\(0.6)(14.7) 


From Table 13 on page 77, selecta 20 inch NPS schedule 40 pipe with d = 18.814 inches 
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© Knowing that D = d/12 = 18.814/12 = 1.568 feet, calculate the relative roughness e /D to be 
9.57 x 10°. Sample Calculation: the Reynolds number for pipe segment AB is: 


W | 370,000 


Ne= 6.31— = 6.31 
* d (18.814)(0.0111) 


| =1,13%10' 
u 


Sample Calculation: knowing these values of ¢ /D and Nre, use a Moody chart, internet applet, or 
any accepted empirical formula to compute a Colebrook friction factor of f= 0.01206. 


© Compute the Mach number at point A using Equation [40] to be: 


M, = 000244951 W, |ZT,  (0.00244951(370,000) |(1(645) _ 0355 


Pd’ VM, (14.7)(18.814)° 55.8 


@ Substitute the calculated and known values for pipe segment AB into Equation [43] to develop 
an expression which can be used to find the value of P;: 


P? 7L If 
l.l m rt n Ta +1 
P d P, 
2 2 
5 —-— U2 100120609) c. " 
P; 18.814 P; 


2 2 
P| = 0.3507|8.577+ n| 2| | 41 
P, P, 


Since this expression is implicit in Ps/P., use iteration to solve for Ps. Begin the iteration process by 
arbitrarily assuming that the pressure at point B is twice that at point A: 


P, = 2P, > P,/P, = 2 


Substitute this trial value for pressure ratio into the test equation and get, 


(2)? = 0.3507 [8.577 im y ] +1 
4 < 4.494 
Try P;/P, = 2.5 to compare the equality: 
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(2.5)? = 0.3507 [8.577 «In (2.5) ] +1 
6250 > 4.651 


Because the inequality sign changed direction, it is obvious that 2 < Ps/P4 < 2.5. 
Try P;/P, = 2.25 to observe the directional change, if any, of the inequality sign, 


5.063 > 4.577 


No change indicates a lower value for P;/P4. See the table that follows for iterative cycles: 


TABLE 11 Iterative Cycles for values of P,/P, to Test the 
Inequality of Equation [43] for Example 12 


Left side Right side 
Pal P4 of Equation Sign of Equation 
2,125 4.516 < 4.537 
2.188 4.787 > 4.557 
2.156 4.648 > 4.547 
2.141 4.584 > 4.542 
2.133 4.550 > 4.539 
2.130 4.537 z 4.538 


As one might imagine, the foregoing exercise could easily be automated somewhat by use of a 
programmable calculator or spreadsheet: 
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E3 Whitesides - Iteration Example 


tB] Ele Edit View Insert Format Tools Data Window Help 
: Calibri -10-|BZ UJ|£E zi $ % 


= (LN(B6)+ 8.577)+0.3507-+1 


The inlet pressure to pipe segment AB is therefore, 


Ps = (2.130)(14.7) = 31.3 Ib/in? absolute 


@ Proceeding upstream and in an identical fashion, a pipe size will now be selected for pipe 
segment BC with P5 now becoming the downstream pressure value. 


0.5 0. 
"DÀ ZT/M,, _ | 185,000 /(1)(598)/46 | _ T 
408.245M P, (408.245\(0.6)\(31.3) 


From Table 13 on page 77, selecta 10 inch NPS schedule 40 pipe with d = 10.020 inches 
For this size, 


Nre = 1.18 x 107 


é/D = 1.79 x 10* 
f = 0.01354 
Mpg = 0.520 


Pc = 38.8 lb/in? absolute 


@ Proceeding upstream and in an identical fashion, a pipe size will now be selected for pipe 
segment CH with Pc now becoming the downstream pressure value. 
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d — 7.981 [8 inch NPS schedule 40] 
Nre = 8.69 x 106 
&gD = 2.26 x 10* 
f = 0.01419 
Mc = 0.426 
Py = 59.6 lb/in? absolute 


Since there is a large positive difference between the MABP and the calculated pressure at point H 
(124.7 versus 59.6), try the next smaller nominal pipe size for pipe segment CH: 


d = 6.065 [6 inch NPS schedule 40] 
Nre = 1.14x 10’ 
&D = 2.96 x 10% 


f = 0.01497 
Mc = 0.737 
Py = 102.9 Ib/in? absolute 


Py < MABP, therefore use the 6 inch NPS schedule 40 pipe size for pipe segment CH. 


See Figure 14b and Table 12 for the remainder of the pressure relief manifold's pipe size 
selections. 


TABLE 12 Summary of Final Data for Pressure Relief Manifold 
Pipe Size Selection - Example 12 


Line AB BD D—F DE B—C C—H CG 
NPS 18 12 8 8 10 6 6 
Schedule 40 40 40 40 40 40 40 
d 18.814 12.000 7.981 7.981 10.020 6.065 6.065 
L JUN 225 17 125 100 300 150 
W 370,000 185,000 60,000 125,000 185,000 110,000 75,000 
T 645 692 800 640 598 610 580 
M, 55.8 71 54 84 46 40 60 
u 0.0111 0.0117 0.0130 0.0110 0.0099 0.0100 0.0098 
MABP  -- --- 45.9 45.7 --- 124.7 58.7 
P; 14.7 31.3 35.7 9597 31.3 38.8 38.8 


M; 0.578 0.320 0.249 0.371 0.520 0.737 0.320 
Nre 1.12 x 10’ 8.31 x 10° 3.65 x 10° 8.98 x 10° 1.18 x 10” 1.14x 10’ 7.96 x 106 
f 0.01280 0.01312 0.01434 0.01418 0.01354 0.01497 0.01500 
P, 3155 Ba 40.2 42.7 38.8 102.9 47.8 


Calculated values are indicated in black; given values are indicated in red. 


Selecting the Optimum Pipe Size 
© 2008, 2015 Randall W. Whitesides, CPE, PE 


www.PDHcenter.com PDH Course M270 


W - 370,000 Ib/hr 


T=645°R 
My = 55.8 
p= 0.0111 cP 
W = 185,000 Ib/hr 
: T=598°R 
Figure M, = 46 


u = 0.0099 cP 


CO;* H;O l 


www.PDHonline.org 


CHi 
W = 110,000 Ib;/hr W = 125,000 Ib/hr 
T=610°R , d- 12 T=640°R 
M; = 40 M, = 84 
Lim 0.0100 cP d= 10 u = 0.0110 cP 
TO © .! oí 
Mo i (c) " d-8 | ERI E & 
a = as 
-2 | dk | K “ye és 
a ^ 4-6 | Eu 
83 TI E 
zs $= 
f d-6 oz 
W = 185,000 Ib/hr 
T2692» R 
w7 1 
g u= 0.0117 cP 
aa 
it (S) 
eo 
rj $ jm — d=8 
ga I SO,* H:0 
=a W = 75,000 Ib/hr 
T =580°R A C.H; 
My = 60 W = 60,000 Ib/hr 
= 0.0098 cP T — T = 800°R 
: d. M, = 54 
(F) u = 0.0130 cP 
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Summary Diagram for Example 12 
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MABP = 45.9 psia 
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Steam Pipes 


Water is not only one of the most common substances and indispensable to life, it is also one of the 
most important media in engineering applications. Steam, the hot gas or vapor phase of water, 
provides process heating, pressure control, mechanical drive, component separation, and is the 
source of hot water for many process reactions." Steam systems can be broken down into the four 
distinct subsystems of generation, distribution, end use, and recovery. Piping transports steam 
between the subsystems. Optimum pipe size determination in the steam distribution and recovery 
components can be particularly challenging in that they offer a unique set of problems. The fluid 
contents of recovery (condensate) piping is a dynamic creature that can have a mixed phase flow of 
a non-Newtonian compressible fluid. Distribution and condensate systems will be dealt with 
separately within this course. 

Steam distribution piping is generally sized with consideration for pressure drop and fluid 
velocity. One simplified approach often used is to arbitrarily select a pipe size and then estimate the 
resulting pressure drop over the length of pipe under analysis, considering steam as an 
incompressible fluid. 


Steam modeled as an incompressible fluid in long lengths of pipe 
can result in inaccurate results because it assumes that the density 


remains reasonably constant. 


The simplified method utilizes a form of the Darcy-Weisbach formula? which has been 
modified to accommodate the insertion of the customary steam flow and specific volume units of 
pounds per hour and cubic feet per pound respectively: 

s f LW?v 


^s [47] 


AP = 336x10. 


Further simplification of Equation [47] with assumed fully turbulent flow, transformed to 
allow for the direct determination of steam pipe size results in, 


d= 175 m [48] 
= Í 2 


Where x is the steam quality expressed as a decimal fraction. 
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Saturated steam is a term used to designate 100% vapor 
saturation. Dry steam is another term for saturated steam. 
Steam quality (x) is the proportion of saturated steam in a 
saturated water/steam mixture (mass fraction of gas phase). A 


steam quality of zero (0) indicates 100% water while a steam 
quality of one (1) indicates 100% steam. Wet steam is a term 
that designates steam with a quality less than one (1). 


Equation [48] is of extreme utility when the recommended steam velocities previously listed 
in Table 8 on page 42 are felt to be applicable. An example follows: 


Example 13 Selection of an Appropriate Steam Pipe Size 
Problem: A relatively short length of 200 Ib/in? steam pipe is to added in a power plant. The line 
will flow 80,000 lb/hr of steam with a moisture content of 9% w/w. Without regard to pressure 


drop, what is the probable nominal pipe size that should be installed to limit noise, vibration, and 
pipe wall erosion? 


Given: P = 200 Ib/in? ; W = 80,000 lb/hr ; wet steam. 


Find: A nominal pipe diameter d for the stated power plant utility application. 


Solution: 

( Without more information, the operating pressure of 200 Ib/in? will be considered high 
pressure steam for the assumption of a suitable fluid velocity. From Table 8, select a median 
flow velocity of 8,000 ft/min. 


@ The stated 9% moisture content translates to a steam quality of X = 1 - 0.09 = 0.91 


Q Referring to any of a number of available steam tables such as Reference 29 (Keenan & 


Keyes), find that the specific volume of saturated steam at 200 Ib/in? gage pressure 
(214.696 Ib/in? absolute) is v= 2.137 ft/lb. 
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® Using Equation [48], the pipe size to limit noise, vibration, and erosion is, 


[Wv x 80,000) 2.137)(0.91 
d= 1754 —— = urs X (0.91) = 7.717 inches 
V 8,000 


®© From Table 13 on page 77 select an 8 inch NPS schedule 40 pipe with d = 7.981 inches 


Once steam has done its work of surrendering latent heat, it condenses to the liquid phase 
referred to as steam condensate. Condensate retains residual useful energy and its physical recovery 
is important for efficient operation. It is returned to the generation unit by means of pumping of the 
hot liquid or through pipe lines at reduced pressures. Pumped condensate return lines in properly 
engineered systems only carry water and as such their pipe size can be determined by the 
conventional incompressible flow methods previously outlined. It is in the low pressure return 
lines, often referred to as headers, where pipe sizing for condensate recovery takes on a new 
dimension: two-phase flow brought about by the larger effective volume and reduced pressure of 
the fluid. When hot condensate under pressure is released to a lower pressure, its temperature drops 
rapidly to the saturation temperature of that lower pressure. The sensible heat released from the 
liquid condensate causes some re-evaporation into flash steam, and the two phases coexist within 
the pipe. Neglecting the presence of flash steam yields undersized pipe sizes that result in increased 
back pressure and velocity, with a liquid phase that could lead to water hammer. Intense water 
hammer can cause damage to personnel, pipes, and equipment. Consequently, the sizing exercise 
must be one which provides some type of two-phase-flow correlation. 

It has been shown that most flow patterns for condensate return pipe lines lie within the 
annular region? although on some occasions they may fall within the dispersed, or even slug-flow 
regions. Generally, the flashed steam constitutes less than 30% w/w of the fluid mixture. Ruskin’! 
developed a method for calculating pressure drop of flashing condensate which is based on the 
assumption of a single homogenous phase of fine liquid droplets dispersed in the flashed vapor. 
This method provides pressure drops that are comparable to those obtained by more rigorous 


methodologies. Using Ruskin's method, Blackwell? 


developed the following flashing correlation 
equations based on the assumption that the vapor-liquid mixture throughout the steam/condensate 
header is represented by the mixture for end-pressure conditions. The subscripts / and g refer to 
liquid and gas, respectively; the subscripts 1 and 2 denote upstream (higher) pressure and 


downstream (low or header) pressure, respectively. See Figure 15. 
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COHDEHSATE RETURN HEADER 


Figure 15 — Relationship between High and Low Pressure Condensate Lines 


The density of the condensate can be correlated to a cubic form in terms of the header pressure: 
p; = - 13x10 P? + A8x10^* P? - 0.078 P, + 60.827 [49] 
The density of the flash steam can be correlated to the header pressure by: 


p, = 0.0029 p^? [50] 


The steam quality of the two-phase mixture is approximated by the expression, 


x = [In A] [0001610621 , 00088] - [o.00671(1n 2^" ] [51] 


Rearranging the original Blackwell? velocity equation allows for the computation of a suitable pipe 


diameter when a mean fluid mixture velocity is known: 
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[52] 
d = 1.7476 


The recommended maximum velocity of flash steam/condensate return lines to limit 
physical pipe deterioration is 5,000 feet per minute with optimum velocities in the range of 3,000 to 
3,600 feet per minute. A worked example which illustrates the use of these equations to determine 
an optimum pipe size for a condensate return header which is expected to have flashing flow 
follows: 


Example 14 Selection of an Appropriate Steam Condensate Return Pipe Size 


Problem: A low pressure steam condensate return header which is designed to operate at 10 Ib/in? 
will receive 10,000 lb/hr of 100 Ib/in? condensate from a process equipment steam trap arrangement 
similar to that shown in Figure 15. Flash steam generation in the header from sensible heat re- 
evaporation is expected. Approximate the probable header pipe size that should be installed based 
on normally recommended fluid velocities. 


Given: P, = 100 Ib/in? ; P; = 10 Ib/in? ; W = 10,000 lb/hr ; flashing conditions. 
Find: A nominal pipe diameter d for the return header portion of the system. 


Solution: 
© Convert the pressures to absolute values: 


P, =100+ 14.7 = 114.7 Ib/in? absolute 
P, =10 + 14.7 = 24.7 lb/in? absolute 


@ Estimate the density of the liquid phase in the header using Equations [49], 


p, = - 13x10 P? - 48x10 * P? - 0.078 P, + 60.827 
p; = - 13x10 5(247)? + 48x107^(24.7) - 0.078(24.7) + 60.827 
p, = 59.17 lb/ft? 
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Q Estimate the density of the vapor phase in the header using Equations [50], 


p, = 0.0029 P?** = (0.0029)(24.7) * = 0.0587 Ib/f? 
€) Approximate the steam quality of the two-phase fluid in the header using Equation [51], 


ee [in B] [o.0001 2822-49 919702) 2 0.0088]- [oo067(n 2*7] 


Il 


[in 114.7} [0.00016120691924 + oooss]- [0.00671(1n 24.7)” | 
0.1044 


X 


Il 


X 


® Calculate the mass flow rate of each phase in the two-phase fluid: 


W - 10,000 Ib/hr 
W, = xW = (0.1044)(10,000) = 1,044 Ib/hr 


W, = W-W, = 10,000- 1,044 = 8,956 Ib/hr 


© Select a suitable fluid velocity based on the recommended range of 3,000 to 3,600 ft/min stated 
in the course: — V = 3,300 ft/min. 


@ Determine a optimum pipe size from Equation [52] based on the calculated data, 


i^ iA 1044. ur 89567, 
Pg JA a 0.0587 59.17 
5 - 17476 


3,300 


d = 1.7476 


= 4.074 inches 


From Table 13 on page 77 selecta 4 inch NPS schedule 40 pipe with d = 4.026 inches 
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Special Pipe Size Design Criteria 


Various industry-specific, non-mandatory special pipe sizing criteria, are often adopted by local 
building codes and thus become regulatory. A good example of a widely adopted method is that of 
the Hazen-Williams formula in the fire protection industry. It is used to calculate pipe size based on 
pressure differential, length of pipe, and pipe material roughness. 


452195 |? : 
div P 
Where Q is in gallons per minute 
The Hazen-Williams formula, written in its original form [NFPA 13 §6-4.2.1], 
4529? 


is more commonly used for pressure drop calculations in sprinkler systems, water distribution 
systems, and irrigation systems. The variable C is a friction loss coefficient which is dimensionless - 
the higher the C factor, the smoother the pipe. Typical C values used in design, which take into 
account some increase in roughness as pipe ages, range from concrete and cast iron at C —100 to 
C= 150 for PVC and FRP. 


While similar in value, the variable C of Equation [53] should not be confused with the 
variable c in the API RP 14E formula, Equation [25]. 


The Hazen-Williams formula is empirical and lacks physical basis. The formula's 
application is severely limited to normal conditions of water flowing at moderate temperatures 
(40° F < T< 75? F) with a mean fluid velocity V ~ 3 ft/sec, Nre > 10°. Even the use of the formula 
with hot water which exhibits a lower viscosity, can result in significant error. The Darcy-Weisbach 
method should be used for other liquids or gases. 
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TABLE 13 Inside Pipe Diameter of NPS Pipe 
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d 
13.624 
13.230 
12.814 
15.500 
15.250 
14.688 
17.500 
17.124 
16.500 
19.500 
19.000 
18.376 
21.500 
21.230 
21.000 
23.500 
23.250 
22.626 
25.250 
25.000 
24.500 
21.230 
27.000 
26.500 
29.250 
29.000 
28.750 
31.250 
31.000 
30.500 
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This course is best summarized by the block diagram shown below: 


OPTIMUM PIPE SIZE 


SELECTION PROCESS 


ENERGY COSTS MOT APPLICABLE 


PRESSURIZED GRAVITY 
FLOW FLOW 


FULL PIPE 


LIFE CYCLE COST (SINGLE PHASE) 


ANALYSIS 


(TWO PHASE) 


ECONOMIC 
ANALYSIS NOT SELF-VENTING 
WARRANTED 


O¥VERRIDING PHYSICALLY OPERATION 
PROCESS SMALL IN TIME 


CONSIDERATIONS SCALE LIMITED 
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